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ABSTRACT 


This  memorandum  records  the  proceedings  of  a  two-part  seminar  series 
sponsored  by  the  Naval  Undersea  Warfare  Center  (NUWC)  Division,  Newport,  Rl 
during  1994.  The  first  part  of  the  seminar  series  dealt  with  the  hydrodynamics  and 
aerodynamics  of  swimming  and  flying  in  nature ;  the  second  part  dealt  with  several 
aspects  of  theoretical  turbulence. 
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FOREWORD 


The  Naval  Undersea  Warfare  Center  (NUWC)  Division,  Newport,  Rl, 
sponsored  a  two-part  seminar  series  in  1994  entitled  "Swimming  and  Flying  in 
Nature."  The  theme  originated  from  a  desire  to  find  out  what  can  be  learned  from 
the  hydrodynamics  of  the  biological  world  that  can  be  applied  to  engineering 
vehicles.  This  series  follows  one  held  in  1992,  entitled  "Turbulence  and  Its 
Control",  and  another  in  1993,  titled  "Microfabrication  and  Turbulence."  The 
proceedings  of  these  preceding  series  appear,  respectively,  in  the  NUWC-NPT  TM 
922089  and  932089 

The  first  part  of  the  1 994  series,  Swimming  and  Flying  in  Nature, 
comprised  four  seminars  devoted  to  the  hydrodynamics  of  swimming  and  the 
aerodynamics  of  flying  in  nature;  the  second  part  comprised  two  seminars  that 
dealt  with  several  theoretical  aspects  of  turbulence.  The  seminars  held  in 
January,  February,  March,  April  and  May  were  an  activity  of  NUWC  Newport's 
Hydrodynamics  Sphere  of  Excellence,  which  is  one  of  NUWC's  leadership  areas. 
The  presentation  materials  used  during  the  seminars,  consisting  mainly  of 
viewgraphs,  are  reproduced  in  this  report  in  their  original  form,  with  only  minor 
reformatting  having  been  done. 

In  the  first  seminar,  Professor  Kashin  gave  the  highlights  of  his  lifelong 
experimental  research  on  the  biological  aspects  of  fish  locomotion,  conducted  in 
Russia  and  Sweden.  He  made  many  interesting  connections  between  the 
physiology  of  a  fish  and  its  hydrodynamics.  Most  importantly,  Professor  Kashin 
examined  the  neural  origin  for  the  generation  of  locomotion.  He  closed  his 
presentation  with  the  videophotograph  of  a  transected  eel,  the  lower  part  of  whose 
body  was  performing  continuous  swimming  movements  when  the  interneural 
network  was  appropriately  excited. 

In  the  second  seminar,  Professor  Newman  gave  a  glimpse  of  his  long 
association  with  the  aerodynamics  of  flight  in  nature,  particularly  of  the  flights  of 
dragonfly  and  nocturnal  owl.  He  brought  attention  to  the  many  features  of  the 
leading  edges  and  of  the  topology  of  the  wing  surface  of  these  flying  animals, 
which  have  no  counterpart  in  engineering  aerodynamics.  Clearly,  animals  in  flight 
adapt  purposefully  to  their  environment  in  ways  that  defy  conventional  engineering. 
It  follows  from  Professor  Newman's  experimental  studies  that  there  are  many  ill- 
understood  mysteries  in  nature  that  are  of  interest  to  hydrodynamics  and 
acoustics. 


v 


Professor  Wu,  in  the  third  seminar,  visited  the  so  called  Gray's  Paradox  on 
the  energetics  of  fish  locomotion.  After  equating  thrust  to  drag,  which  is  mostly  the 
pressure  drag  for  a  fish,  Professor  Wu  used  potential  flow  analysis  to  conclude 
that  the  drag  of  a  live  undulating  fish  is  significantly  lower  than  that  of  a  freshly 
dead  fish,  or  the  drag  of  a  rigid  fish  model.  Professor  Wu  proposed  that  the  origin 
of  this  economical  locomotion  in  a  fish  lies  in  unsteady  hydrodynamics.  It  is 
tempting  to  think  that  the  study  of  the  unsteady  hydrodynamics  of  aquatic  animals 
might  lead  to  the  development  of  new  concepts  on  the  locomotion  of  engineering 
vehicles. 

The  fourth  seminar  by  Professor  Lumley,  unfortunately,  was  canceled  due  to 
the  sudden  illness  of  his  wife. 

The  second  part  of  the  seminar  series  was  brief,  consisting  of  only  two 
lectures.  In  the  first,  Professor  Walker  dealt  with  the  transport  of  turbulence  vorticity 
away  from  a  solid  surface.  He  used  inviscid  vortex  analysis  to  argue  the  existence 
of  very  thin  shear  layers  that  are  normal  to  the  wall  in  a  turbulent  boundary  layer. 
Professor  Walker's  work  suggests  that,  in  a  time-dependent  viscous  computation, 
very  fine  grid  resolution  is  needed  to  capture  any  intermittent  outward  eruptions  of 
vorticity  from  the  wall. 

In  the  second  seminar  on  theoretical  turbulence,  Professor  Sreenivasan 
showed  how  one  can  build  the  synthetic  time  series  of  any  turbulence  variable  that 
has  the  "right"  statistical  properties.  He  has  distilled  the  knowledge  of  the  fractal 
nature  of  turbulence  to  invent  this  method  of  creating  a  synthetic  time  series.  The 
utility  of  his  invention  seems  to  lie  in  the  time  savings  that  might  be  derived  from  its 
use  as  the  initial  condition  for  generating  a  turbulence  data  base  using  the  so- 
called  Direct  Numerical  Simulation  method. 
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Swimming  and  Flying  in  Nature 


1.1  Fish  Swimming 


Sergei  Kashin 
Northeastern  University 
(Formerly  of  Moscow  State  University) 
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NUWC  Division  Newport 

SEMINAR  NOTICE 


FISH  SWIMMING: 

KINEMATICS,  NEURAL  AND  MUSCULAR  CONTROL 

Dr.  Sergei  Kashin 

Department  of  Biology  and  Marine  Science  Center 
Northeastern  University 
(formerly  with  Moscow  State  University) 


20  000  fish  species  show  great  diversity  in  modes  of  swimming. 
The  most  common  mode  is  the  undulation  of  the  body  or/and  fins.  The  wave 
of  undulations  travels  along  the  body  and  usually  low  speed  is  controlled 
by  the  amplitudes  of  undulations  whereas  high  speed  is  controlled  by 
frequency  of  undulations.  The  fish  musculature  is  divided  onto  a  number  of 
segments  along  the  body.  The  number  of  segments  is  exactly  equal  to  the 
number  of  vertebra  and  each  segment  or  myotome  is  supplied  with  its  own 
pair  of  sensory  (dorsal)  and  motor  ventral)  nerve  roots  from  spinal  cord. 
The  structural  and  functional  organization  to  decrease  the  number  of 
degrees  of  freedom  of  system  is  discussed.  The  nervous  commands  to  the 
muscular  system  were  investigated  in  fish.  The  functioning  of  “locomotor 
center”  in  fish  brain  is  described  as  well  as  its  relation  to  the  spinal 
generators  of  rhythm.  The  behavior  of  molecular  units  of  contraction 
(sarcomeres)  was  studied  in  intact  fish  with  laser. 

Some  applications  to  the  understanding  the  hunting  behavior  as  well 
as  some  problems  in  engineering  a  fish-like  robot  are  discussed. 


Thursday,  the  27th  January  1994 
Conference  Room,  Bldg.  679  First  Floor 
Time:  10:30  AM 
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T0PI1CS  TODAY 

Movements  of  the  body  to  produce  swimming 

Muscle  architecture 

How  muscle  work  to  produce  movements 
Nervous  control  of  movements 
Molecular  contractile  activity  in  vivo 
Fish-like  robots 


1.1-5 


1.  FORM,  FUNCTION,  AND  LOCOMOTORY  HABITS 


2 


panded,  is  followed  here.  Examples  of  fish  displaying  these  modes  are 
shown  in  Fig.  1.  As  Breder  stated,  the  suffix  “-form”  (e.g.,  in  anguil- 
liform)  refers  to  the  types  of  movement  and  not  to  the  body  forms,  and 
is  therefore  not  strictly  parallel  to  words  such  as  “fusiform.”  Indeed, 
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Fig.  I.  Mode*  of  forward  sw  imming  in  fish,  arranged  alone  the  vertical  axis  according 
to  the  propuhise  contributions  of  !>od>  and  fins  indicated  b>  density  of  shading),  and 
alone  the  hnri2ont.il  axis  according  to  a  scale  running  from  serpentine  undulation  (more 
than  one  wavelength  »•' 
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a  -  movement  of  aniomal 
forward 

b  -  movement  of  bending 
backward 
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CARANGIFORM  MODE 
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The  evolution  of  locomotor  system  of  fish  could  be  considered  as  a  process 

of  multipliction  of  a  number  of  joints  (the  degrees  of  freedom)  or  ... 
as  reduction  of  it. 
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Myotome  ( myomere )  -  al  I  musle  fibers  between 
two  neghboring  myosepts.  Each  myotome  has  its  own 
innervation  from  spinal  cord  -  two  dorsal  and  two 
ventral  spinal  nerves. 
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Backbone  Structure 


Dimensions  (I).  Scaled  to  the 
bodylength  (L).For  mackerel  or 
perch  type  l=0.66L/N,  where  N  is 
number  of  vertebras.  And  1=  L/N  for 
long  fish  like  lamprey  or  eel. 

Number.  Equal  to: 

-  24-30  for  bass  and  perch 

-  30-40  for  mackerel  type  (tuna-39) 

-  40-60  for  cod,  herring,  salmon  etc. 

-  >100  for  lamprey  and  eel 

Shape.  Short  cylinders  (d/l=1)  with 
attached  processes  and  arcs.  In 
vertical  plane  -  neural  and  hemal 
arches  and  processes.  Processes 
are  used  as  levers  for  rotarion. 
Horisontal  -  ribs  -  are  also  used  for 
muscle  attachment.  Zygapophysis 
processes  restricting  rotation. 

Material.  Bone  consists  of 
polysaccharide  and  protein  matrices 
where  calcium  phosphate  is 
deposited. 


DORSAL 
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VENTRAL 
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Angles  of  muscle  fibers  vitn  the  main  axes  of  the  body  as  a  function  of 
the  cross-sectional  co-ordinates  (upper  right). 

Projections  on  the  horizontal  plane  (upper  left) 

Projection  on  the  vertical  plane  (belov) 
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Pmc.  4.  Ilojiamw*  toh®k  cptAmft  jutxxx  ma  b  nocjt«AOBJtT«JKM&a 
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Lateral  displacement  of  different  parts  of  the  body  of  the 
eel. 
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Speed  of  swimmind  vs  frequency  and  periodlength 
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Frequency  of  undulations  vs  speed  of  swimming 
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The  tailfin  amplsvvimming  speed 
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Three  types  of  the  distribution  od  amplitudes  along  the 
fish  body  (1  -  eel,  2  -  trout,  3  -  sculpin) 
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Vb 


Velocity  of  bending  wave  along  the  body  vs  speed  of 
swimming  for  three  species  (goldfish,  eel  and  trout) 
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Electromyography  of  the  Ciprinus  carpio 
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Electromyography  of  trout 
Midline  -  white  muscles 
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for  three  speeds  of  swimming. 


* 


A  0.15  M/« 


B  0.41  m/c 
0  500  M/c 


O  t  « 


C  063  M/c 

0  250  MC  0  10 


D  117  M/C 

0  Z50 MC  O  i  ® 


F <5r  graphs  represents  the  EMGs  recorded  from  eel  swimming  wltK 
four  different  speeds. 

Bars  on  the  left  side  of  each  graph  represent  the  duration  of  EMG 
(Dots  show  the  cycle  duration). 

The  faster  fish  swims  the  shorter  the  delays  along  the  bodylength. 
However  the  phase  is  constant  as  it  is  seen  on  the  right  side  o^ 
each  graph. 
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Stimulation  of  the  locomotor  center  in  the  midbrain. 
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Thrust  Locus 


1.1-28 


Fig.  2.  Propagation  of  EMG  waves  along  the  body  of  lamprey.  Based  on  Williams  et 
1989.  All  temporal  and  linear  parameters  are  normalized  to  the  bodylength 

(L)  and  the  periodlength  (T). 
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MARINE  SYSTEMS  ENGINEERING  LABORATORY 


UNDULATING  ROBOT  PROGRAM 


BIOLOGICAL  RELATIONSHIPS  IMPORTANT  TO 
ENGINEERING  MODEL  DESIGN 


•  Time  delay  between  activation  of  segments  is  always 
a  fixed  percentage  of  cycle  time. 


•  Frequency  of  oscillation  is  linearly  related  to  speed. 


4ft  Increasing  stiffness  is  required  to  increase  frequency 
hence  speed. 

•  Left  and  right  sides  of  system  are  alternately  activated. 


•  Duration  of  excitation  does  not  exceed  1/2  of  cycle  time. 


m  Different  muscles  are  used  for  slow  and  fast  speed 


m  Thrust  is  in  form  of  pulses  at  twice  the  frequency. 


•  Amplitudes  of  undulation  are  directly  correlated  to  basic 
structure,  ie  Mass,  cross  section  area,  and  areas  of 
lateral  projection. 


NORTHEASTERN  UNIVERSITY 


1.1-30 


f  MARINE  SYSTEMS  ENGINEERING  LABORATORY 

I  UNDULATING  ROBOT  PROGRAM 

If  ROBOT  VS  FISH  PROPULSION 


ROBOT: 

Local  flexures  propagate  from  front  to  rear  pushing 
water  back  and  thus  produces  thrust. 

Movement  of  tail  from  side  to  side  also  produces  thrust 

It  will  have  2  basic  frequencies:  1  hz  for  slow  swimming 
and  10  hz  for  fast  swimming. 


»  FISH: 

Bending  waves  propagate  along  body  from  front  to  rear 
generating  thrust.  In  good  swimmers  fish  moves  about 
0.8  bodylength  in  1  cycle. 

The  frequency  of  undulation  varies  smoothly  from 

0.5  to  20  hz  (25  cm  fish).  Larger  fish  have  lower  maximum 

frequency. 

The  amplitude  of  the  trailing  edge  of  the  tail  is  maximal 
and  equal  to  0.2  bodylengths  for  all  species  and  almost 
all  ranges  of  speed. 


NORTHEASTERN  UNIVERSITY  r 
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MARINE  SYSTEMS  ENGINEERING  LABORATORY 


BIOLOGICALLY  BASED  UNDULATORY 
AUTONOMOUS  UNDERWATER  VEHICLE 


FISH  VERSUS  ROBOT  MUSCLE  ANALOGY 


r 

Parameter: 

Fish 

Nitinol  1 

Muscle  Contraction  .  .  . 

.  .  .2% 

3%  -  5% 

Frequency  . 

.  .  .  0.5  -  20hz 

0.1  -  10  hz 

Nominal  efficiency  .  .  . 

....  30% 

5% 

Muscle  Types . 

.  .  .  slow/fast 

slow/fast 
variable  sizes 

Force-velocity  rel . 

...  inverse 

inverse 

Max.  Cycles . 

.  .  .  lifetime 

several  million 

.  _  _ J 

SOME  MAJOR  ENGINEERING  ISSUES 

*  Optimize  number  of  segments  to  achieve  smooth  motion 


*  Minimize  size  and  number  of  SMA  wires  to  minimize  energy  needs 

*  Control  and  coordination  of  multiple  elements:  maneuverability, 
rigidity,  speed,  and  buoyancy 

*  Control  of  temperature  (heating  &  cooling)  to  maximize  frequency 
response 

*  Mechanical  design  &  fabrication  of  skeleton,  mass  distribution 
system,  &  mechanism  to  transfer  internal  motive  forces  to 

interface  with  water  column 


NORTHEASTERN  UNIVERSITY  i 


-  Complicated  muscular  morphology 
results  from  the  necessity  to  bend  the  body  in 
any  particular  point 

-  Red  muscle  system  functions  in  hydro¬ 
static  skeleton,  whereas  white  muscle  belongs 
to  the  lever  type 

-  Cinematically  low  speed  is  controlled  by 
the  tail  amplitude  and  high  speed  by  the 
frequency  of  undulations 

-  Basic  synergies  are  constant  in  phase. 
Only  period  length  is  controlled  by  the  brain 
locomotor  center.  Recruitment  (stiffness  of 
the  body?)  and  intersegmental  phase  lag  are 
controlled  by  the  frequencies  perhaps  on 
spinal  level 

-  All  that  simplifies  the  problem  of 
designing  the  fish-like  robot 


1.1-33/1.1-34 
Reverse  Blank 


1.2  Aerodynamics  of  Bird  Flight 


Barry  G.  Newman 
McGill  University 
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Newport  Seminar  Notice 


LOW-SPEED  AERODYNAMICS:  ANIMAL  FLIGHT  AND  SKIN  FRICTION 

Professor  B.  G.  Newman 

Canadair  Professor  of  Aerodynamics 
McGill  University,  Montreal,  Canada 

As  suggested  by  Dr.  Bandyopadhyay,  three  studies  will  be  discussed: 

1.  Large  Dragonflies  have  a  wingspan  of  90  mm  and  an  unusual  wing  section  consisting 
of  pleats  near  the  leading  edge.  Hand-launched  gliders  with  this  section  were  tested 
at  Re  a  5  x  104  and  showed  that  the  maximum  lift  to  drag  ratio  was  about  10, 
slightly  better  than  that  for  arc  aerofoils  often  used  on  model  aircraft  at  this  Reynolds 
number.  The  separation  and  reattachment  of  the  flow  in  the  pleats  generates  a 
turbulent  boundary  layer  which  increases  the  maximum  lift. 

Other  aspects  of  wing  geometry  will  be  mentioned. 

2.  Nocturnal  owls  have  wings  with  at  least  three  unusual  features  that  are  thought  to 
contribute  to  their  quiet  flight.  Two  of  these  are  leading-edge  combs  and  downy 
surfaces. 

The  combs  are  formed  by  a  curved  extension  of  the  barbs  on  the  exposed  leading  edges 
of  primary  feathers.  They  probably  produce  trailing  vortices  which  inhibit  separation 
near  the  leading  edge  (the  common  Alula  also  does  this  but  on  a  larger  scale).  Smoke 
tunnel  studies  show  that  a  single  primary  feather  has  separation  suppressed  for  about 
6°  incidence  due  to  the  presence  of  the  comb  at  Re  «  104.  Suppressing  the  high 
pressure  fluctuations  near  reattachment  reduces  noise  propagation. 

Down  (1  mm  long  hairs,  10  p  diameter)  is  found  on  the  top  of  the  wing.  Flat  plate 
boundary  layer  measurements  behind  downy  surfaces  showed  that  the  skin  friction 
coefficient  Cf  was  increased  by  about  85%  :  however  the  downy  surface  was  rough 
but  never  fully  rough.  In  contrast  the  longitudinal  turbulence  was  reduced  very  close 
to  the  surface  and  with  it  the  surface  pressure  fluctuations.  Thus  down  is  an  unusual 
type  of  roughness,  not  bluff  and  turbulent,  but  streamlined  and  laminar. 

3.  An  instrument  has  been  developed  for  measuring  skin  friction  on  a  “rough”  surface. 

dU 

Three  pitot  tubes  are  used  to  measure  the  velocity  gradient  —  at  two  places  within 

ay 

the  wall-low  region  where - r  — —  =  —  a  constant,  independent  of  viscosity 

c-)’  dy  K 

and  the  geometry  of  the  roughness.  The  origin  for  the  wall  distance  y  is  uncertain, 
which  explains  why  three  pitot  tubes  are  required. 

The  final  formula  for  the  skin  friction  tw  is  precise  and  the  instrument  does  not 
require  any  calibration.  Measurements  have  been  compared  with  a  Dickinson  skin 
friction  balance. 

Tt  could  be  used  for  ribietted  surfaces  as  well  as  those  with  marine  growth. 

Thursday,  the  10th  February  1994 
Conference  Room,  Newport  Bldg.  679  First  Floor 
Time:  10:30  AM 
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NUWC  Seminar 


Feb.  10,  1994 


Newport,  R.I. 


Dr.  Bandy opadhyay  has  asked  me  to  review  some  past  and  present  work  on  animal 
flight  and  the  measurement  of  skin  friction. 

k 

I  will  first  talk  about  dragonfying,  the  aerodynamics  of  anizoptera-odonata,  and 
then  describe  a  new  robust  way  of  measuring  skin  friction,  which  leads  into  the  final 
topic,  the  aerodynamics  of  nocturnal  owls.  Presumably  some  of  the  "devices"  which 
have  worked  (i.e.  been  developed  on  these  animals)  in  air  might  be  worth  trying  in  water. 

The  work  I  am  going  to  describe  is  mainly  curiosity-motivated. 


1.  Dragonfly  Hertel  Structure,  Form  &  Movement  (1966) 

2  Slides  -*  1  overhead  Pleated  aerofoil 

Structure  and  model  for  smoke  tunnel 

R  =  io4  very  difficult  to  generate  turbulent  b.l. 
to  keep  the  flow  attached. 

2  Model  gliders  were  tested 

Re 

/  Large  4  x  104 
(Small  1  x  104 


plus  McBride  B7 

(40  chords  travelled  for  every  flap 
:  25  Hz  is  frequency 
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27.  Model  Tests  on  a  Wing  Section  of 
an  Aeschna  Dragonfly. 

B  G  NEWMAN 

Department  ol  Mechanical  Engineering  McGill  University.  Montreal  Canaria 

S  B  SAVAGE  AND  D  SCHOUEllA 

faculty  ot  engineering  McG'H  University  Montreal  Canada 

ABSTRACT 

[he  chord  Reynolds  number  ol  a  dragonfly  wing  Hying  al  high  \peed  i\  ol  Ihe 
order  of  Id4  It  is  well  known  that  aerofoils  at  moderate  incidence  experience  an 
increase  ol  drag  coefficient  and  a  decrease  ol  lift  coefficient  as  the  Rcsnolds  number 
is  reduced  below  a  critical  value  of  about  *  It)4  and  that  this  is  associated  with 
complete  separation  ol  the  laminar  boundary  laser 

I  he  aeroloil  section  ol  the  dragonfly  is  unusual,  as  indicated  by  I  ig  la  which  is  a 
skcti  h  of  a  photograph  ol  the  section  gist  inboard  ol  the  nodus  on  the  w  ing  ol  4  cm  four 
/.  r, -1111111 

Attention  is  drawn  to  the  minute  saw-teeth  on  each  web  ol  the  7  section  which 
forms  the  leading  edge  (costa),  and  the  spurs  on  both  sides  ol  the  matrix  which 
supports  the  rear  membrane.  The  saw-teeth  very  likely  act  as  turbulators  to  promote 
transition  in  the  separated  shear  layer  and  subsequent  reattachment  of  the  boundary 
layer  in  a  turbulent  state  Smoke  tunnel  studies  of  the  wing  section  support  this 
hypothesis.  Trapped  vortices  arc  observed  in  the  k'  sections  at  low  incidence,  and  at 
high  incidence  labout  10  I  the  flow  is  observed  to  reattach  to  the  rear  cambered 
membrane,  the  leading-edge  separation  bubble  is  then  greater  than  half  a  chord  in 
length.  The  purpose  of  the  spurs  which  are  roughly  the  height  of  the  viscous  sub-layer 
for  the  reattached  flow,  is  less  obvious 

The  average  flapping  frequency  in  forward  flight  labour  10  m  s  1 1  is  roughly  .  k  Hr 
(Nachngall" Jl)  so  that  in  one  cycle  the  wing  moves  forward  about  40  chord  lengths. 
It  is  therefore  postulated  that  the  aerodynamics  may  be  usefully  studied,  at  least 
initially,  on  a  static  wing  in  steady  flow. 

Indoor  free-llight  tests  have  been  made  on  two  sizes  of  model  glider  with  a  similar 
wing  section.  Under  stroboscopic  illumination  flights  have  been  photographed  using 
a  long  time  exposure.  With  reflectors  placed  on  the  model  and  above  the  horizontal 
floor  the  flight  altitude,  flight-path  angle  and  speed  have  been  measured  in  steady 
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PRACTICAL  APPLICATIONS  OF  TWO-DIMENSIONAL 
MEMBRANE  RESEARCH 


B.G.  Newman 

Professor  of  Aerodynamics,  McGill  University 
Montreal,  Canada 


ABSTRACT 


NOMENCLATURE 


Theories,  mainly  inviseid,  were  developed  several  years  ago 
for  two  dimensional  lifting  membranes.  More  recently  experi¬ 
ments  have  shown  that  the  predictions  may  be  very  inaccurate 
for  medium-to-large  camber  ratios.  Local  separations  of  the 
Sow  are  the  cause  of  these  discrepancies. 

Although  two-dimensional  this  work  has  some  relevance  to 
three-dimensional  sails  tacking  obliquely  into  wind.  For  exam¬ 
ple  it  emphasises  the  importance  of  eliminating  leading-edge 
separation  on  both  head-and  main-sails.  Prediction  of  the 
aerodynamic  forces  on  a  three-dimensional  sail  with  separa¬ 
tion  bubbles  is  very  difficult  and  has  not  yet  been  attempted. 

Sail  efficiency  can  be  improved  by  the  judicious  addition 
of  stiffness,  usually  with  battens.  Membrane  theory  has  been 
extended  to  cases  where  the  stiffness  may  vary  from  leading 
to  trailing  edge,  so  that  design  for  a  particular  pressure  distri¬ 
bution  becomes  possible.  Again  only  flattish  sails  give  results 
which  agree  reasonably  with  the  theory. 

Luffing  instability  is  experienced  on  membranes  as  well  as 
sails  and  is  shown  to  be  a  function  of  incidence  and  camber. 
According  to  linear  theory  it  is  a  divergent  instability  and  oc¬ 
curs  just  after  the  sail  becomes  S  shaped  at  a  slightly  negative 
incidence. 

For  square  sails  and  spinnakers  two-dimensional  membrane 
theory  indicates  how  to  increase  drag  and  reduce  oscillations. 
Flat  sails  are  preferable  although  the  increases  in  tension  may 
be  unacceptable. 


e  length  between 

Ci  lift  coefficient 

CPm  base  pressure  coefficient 

Cj  leading-edge  suction  coefficient 

Ct  tension  coefficient 

E  Young’s  modulus  of  elasticity 

I  ■  second  moment  of  the  cross  section  of 
a  stiffened  membrane 
l  length  of  membrane 

n  frequency  of  oscillation 

q  dynamic  pressure  of  the  free  stream 

Ue 

R,  Reynolds  number  — 

1/ 

S  suction  force  at  the  leading  edge  in 

approximately  the  negative  x  direction. 

U  free  stream  velocity 

x  along  the  chord 

y  perpendicular  to  the  chord 

a  incidence 

7  camber  ratio 

i  i  .  I  —  e 
c  excess  length  ratio,  - 

v  kinematic  viscosity 


membrane  supports 

L 

\plPc 
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BOOMERANGS 


by  Barry  G.  Newman,  FRAeS 
Professor  of  Aerodynamics 
McGill  University,  Montreal 


Dr  Barry  G.  Newman  was  educated  at 
Manchester  Grammar  School  and 
Cambridge  University.  He  completed  his 
PhD  at  Sydney  University  in  1952.  He 
has  been  the  Canadair  Professor  of 
Aerodynamics  at  McGill  University  since 
1959.  Previously  he  was  Lecturer  in 
Aeronautics  at  Cambridge  University  and 
before  that  he  undertook  flight  research 
with  both  the  Royal  Australian  Air  Force 
and  the  National  Aeronautical  Estab¬ 
lishment,  Canada.  He  has  broad 
research  interests  in  aerodynamics 
ranging  from  turbulent  shear  flows, 
wakes,  jets  and  vortices,  boundary  layers 
and  boundary-layer  control,  Coanda 
effect  and  the  aerodynamics  of  animal 
flight,  air-cushion  vehicles,  wind  turbines 
and  large-scale  flexible  membranes 
such  as  sails,  parachutes  and  inflated 
buildings.  He  is  currently  a  consultant  to 
Pratt  &  Whitney  Canada  Ltd  and  the  Pulp 
&  Paper  Research  Institute  of  Canada. 
He  is  a  Fellow  of  the  Canadian 


Aeronautics  and  Space  Institute,  the 
Royal  Aeronautical  Society  and  the  Royal 
Society  of  Canada. 


Returning  boomerangs  are  used  by 
certain  tribes  of  aborigines  in  Australia 
and  appear  to  have  been  invented  by 
them.  In  this  paper  both  aboriginal  and 
“modern"  aerodynamic  boomerangs  are 
described.  The  flight  dynamics  of  return¬ 
ing  boomerangs  is  qualitatively  ex¬ 
plained  with  particular  reference  to  the 
work  of  Walker,  Lanchester  and  Hess. 
The  non-dimensional  parameters  which 
determine  the  main  features  of  each  part 
of  the  motion  are  identified,  and  simplified 
analyses  are  presented.  The  predicted 
behaviour  is  compared  with  boomerang 
flights.  The  predictions  include  the 
connection  between  boomerang  in¬ 
clination  and  advance  ratio  at  launch  and 
the  effect  of  scale  and  density  on  range 
and  on  the  final  hovering  mode. 

1.2-11/1.2-12 
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Photo:  The  Australian  Information  Service. 
London. 


Aboriginal  Boomerangs 

There  are  three  main  types  of  Australian 
boomerang:  the  non-returning  killer  or 
war  boomerang  which  is  usually  asym¬ 
metrical  with  arms  of  different  length,  the 
returning  boomerang  which  is  usually 
lighter,  and  the  ceremonial  boomerang 
which  is  not  thrown.  Non-returning 
boomerangs  are  not  peculiar  to  Australia 
and  they  are  traditionally  found  on  at  least 
three  other  continents.  Sometimes  these 
non-Australian  weapons  are  made  of 
metal  in  the  form  of  throwing  knives  and 
pictures  of  them  may  be  found  in  Ref  i 
Although  some  of  these  weapons  do  have 
a  tendency  to  return,  the  fully  returning 
boomerang  appears  to  have  been  an 
invention  of  the  Australian  aborigine ,2-3 
Two  examples  of  asymmetrical  non¬ 
returning  Australian  boomerangs  are 
shown  among  the  author's  collection,  Nos 
1  and  2  in  Fig  1 .  Sometimes  a  hook  is 
added  to  the  short  end  of  the  boomerang 
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DRAGON  FLY  EXPERIMENT 
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EXPERIMENT 
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DRAGON  FLY  EXPERIMENT 
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DRAGON  FLY  EXPERIMENT 
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DRAGON  FLY  EXPERIMENT 
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DRAGONFLY  WING  SECTION 


Me  BRIDE  B-7 
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0.08 


0.06 


0.04 


0.02 


LARGE  MODEL 

(Re  ~  3.5  xIO4— 6  x!Q4) 


SMALL  MODEL 
(Ke^l.2  x  I04-2.5xI04) 


DRAGONFLY'  SECTION 
(NO  TURBULATORS) 
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Nocturnal  Owls 


as  distinct  from,  say  Snowy  Owls 

Hunt  by  diving  on  rodents  at  night. 

Silent  flight  is  usually  attributed  to 

(a)  Downy  surface  on  TOP  and  biggest  near  trailing  edge  of  the 
downstream  primary  feather. 

(b)  Shawl-like  trailing  edge. 

(c)  Comb  at  exposed  leading  edge  of  the  primaries. 
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FIGURE  2  {right):  Primary 
feather  of  a  Great  Horned 


Flat  plate  measurements 


downy  surfaces  surrounded  by  silk  velvet  (Fig^.  2  a»d  £  i-f/t). 
Pitot  and  hot  wire  traverses. 


U  ,  yUT 

—  /  Jo0io  -  =  5.6  +  -  -  -  Figs.  3-5. 

Ur  V 

UT  ->  Cf 


iUT 

Down  increases  Cf  by  about  85%,  but  not  independent  of - ,  i.e. 


i/ 


not  “fully  rough”. 

O.K.  for  steep  descent. 

Longitudinal  turbulence  (Fig. 6),  downy  surface  has  reduced 
yJJT 

turbulence  for  - — -  <  100.  Similar  to  surface  pressure  fluctuations. 


v 

Down  ‘rough’  due  to  laminar  flow  over  the  down  in  contrast  to 
conventional  roughness. 

Spectra  shows  a  shift  to  higher  frequencies.  10  x  about.  Fig.  7. 
This  could  then  be  attenuated  by  long  grass  and  low  shrubs. 
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Leading-edge  comb. 

The  leading  edge  of  all  exposed  leading  edges  of  primary  feathers 
has  a  comb  which  is  an  extension  of  the  barbs. 

We  compared  leading-separation  bubble  for  such  a  feather  with  flat 
plate  in  a  smoke  tunnel.  Slightly  dependent  on  Reynolds  number 
(based  on  feather  chord).  Fig.  '• 

Comb  on  feather  suppresses  the  bubble  by  A  oc=  6°. 

x  R 

(note  —  increases  only  slightly  with  roughness) 
c 

No  comb  1  rr™  =  0.12  at  reattachment  cf  0.03  no  separation. 

oo 

0.09  at  reattachment  if  downy. 
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Position  of  tlmo  mean  reattachment  zone  (fully  smooth  surface) 


SMOOTH 
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Table  1.  Boundary  layers  on  the  large  flat  plate. 


SURFACE 

U 

00 

m/s 

6X 

mm 

cf 

U 

T 

u 

oo 

AU 

Ut 

k  mm 
s 

Equivalent 

Sand 

Roughness 

k  U 
s  T 

V 

k  mm 
Measured 
Roughness 

13.4 

8.24 

.0035 

.042 

SMOOTH 

27.8 

7.00 

.0030 

.039 

13.3 

9.85 

.0053 

.051 

5.40 

.810 

37 

.57 

ROUGH 

27.1 

9.58 

.0052 

.051 

8.00 

1.120 

75 

.57 

13.7 

.  9.08 

.0059 

.054 

6.00 

.9 

DOWNY 

30.4 

8.52 

.0058 

.054 

8.35 

.9 

SILK  VELVET 

14.0 

9.23 

.0056 

.053 

5.35 

.7 

0  0.5 
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3.  Measurement  of  Skin  Friction 


An  instrument  consisting  of  3  pitot  tubes  mounted  in  the  wall  law  region  of  a 
turbulent  B.L.  can  be  used  to  obtain  the  skin  friction  for  both  smooth  and  non 
smooth  walls.  It  is  a  sort  of  extension  of  Preston’s  method. 


Rotta 

jT  _  V_  9U  1 

r  V  P  '  UT  dy  k 

ind.  of  viscosity  and  all  aspects  of  the  roughness,  geometry,  size,  if  not  too  big. 


Number  pitots 


Ui  +  y-i  (  U\  ~~  U2 
2Ut  V  2/1  —  2/2 


2  2 

and  2  -  3  yi  -  Vi  =  <*12,  Ui  -  U2  =  ,  ■  .  (Pi  -  P2)  =  -7 yf  "Jr  V  Al2 

p[Ui  +  U2)  p[U\  +  U2) 


=  formula  in  terms  of  A12,  A23  &:  a23 


Checked  for  smooth  walls 

and  #40  and  #24  sandpaper  roughness. 

Pitots  project  sufficiently  ahead  of  support,  and  static  tube  is  mounted  to  the 
side.  Compared  with  skin  friction  measured  on  a  floating  element  balance 
Null  -  LVDT  linear  voltage  differential  transducer 
i.e.  variable  reluctance  transducer 


Possible  error  atJ-  O.K.  3D  analysis  using  sources  I  unimportant 
Pressure  difference  A, j  Low  Re  Dynamic  8,  10,  11,  13 

Gaussian  smoothing  Barocel  590  Capacity  type 

Fig.  10 

±  •  IPa  — ►  <  1  Hz  sequential  readings 
Possible  use  Hull  roughness  due  to  marine  growth. 
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Stainless  Steel  Hypodermic  Tubing 
0.82  ±  0.01  mm  dia. 


K 


f —  Brass 
/  1 .59 


Tubing 
mm  dia. 


mm 


ai 2  —  1  -57  ±  0.01  mm 
a2^  =  1 .62  ±  0.01  mm 

Figure  2:  Three  Tube  Instrument 


1.2-52 


1.2-53 


lollowing  Kotta  (1962),  if  the  surface  roughness  is  not  too  large, 


y_dU  =  1 

UT  dy  /c’ 


where  UT  =  (tw/p)i  is  the  skin  friction  velocity  and  k  «  0.41  is  the  von  Karman 
constant.  Approximating  the  velocity  gradient  at  the  midpoint  between  tubes  1 
and  2  by  a  linear  interpolation  of  the  velocity, 


y\  +  y2  /  U\  ~  \  _  1_ 

2UT  V  yi  -  y2  y  K 


provided  yi  —  y2  yi-  Similarly,  for  tubes  2  and  3, 


Vi  +  ys  fUi-  Us\  _  3.  ,  . 

2UT  \y2-yj~  k'  U 

provided  y2  —  y3  <  y2.  The  origin  for  the  y-coordinate  is  difficult  to  establish 
for  a  rough  surface  and  therefore  y  is  eliminated  by  taking  the  difference  of  these 
two  equations.  This  leads  to  the  expression 


J_  _ _ 2 _  /  o-n _ a23  \ 

UT  /c(fll2  +  fl2s)  —  U2  U2  ~  U$) 


where  al2  =  yi  -  y2  and  e23  =  Vz  ~  ys- 

Since  pressures  rather  than  velocities  are  measured,  the  velocities  are  ex¬ 
pressed  as 
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Figure  1:  Three  Pitot  Tubes  in  the  Logarithmic  Region 

Expanding  the  velocities  to  second  order  in  Ay/ (ft  —  p)  and  substituting  into 
Eq.  5  gives 

1  _  [fy>(ft  ~  P)}1/2  0,12  _  <*23  Ol2  +  <*23 

Ur  /c(ai2  +  <123)  .A12  A23  4(^2  —  p) 
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Generalized  slender-body  theory  has  been  fruitfully  applied  to 
the  bio-hydrodynamical  analysis  of  modes  of  locomotion  of  elongated 
fishes  by  means  of  body  undulation.  Of  interest  are  the  highly 
efficient  modes  of  carangiform  (employing  movements  of  only  fish's 
posterior  end),  thunniform  (marked  by  a  lunatail  caudal  fin), 
balistiform  (propagating  distally  waves  along  dorsal  and  anal  fins 
with  body  held  rigid)  and  their  derivatives.  Studies  of  the  vortex 
dynamics  involved  in  these  modes  of  swimming  motion  have 
elucidated  the  basic  mechanism  of  interaction  between  the  vortex 
sheets  and  the  body-fin  system,  pointing  out  how  thrust  can  be 
enhanced,  loss  of  hydromechanical  energy  reduced,  and  biodynamic 
efficiency  improved.  Also  discussed  will  be  the  problem  of 
energetics  of  fish  swimming  from  a  mechanophysiological  approach. 
From  a  series  of  comparative  studies  of  the  scale  effects  of  fish 
locomotion,  Gray's  paradox  appears  to  be  substantially  clarified. 
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Fig.  7.  Swimming  endurance  of  sockeye  salmon  and  rainbow  trout.  Transition  zones  between 
different  levels  of  activity  are  shown  with  broken  lines.  (After  Brett',,,JJ|). 


V,  Swimming  speed  (m  $"') 

Fig.  8.  Specific  energy  cost  S  as  related  to  swimming  speed,  for  adult  sockeye  salmon.  Ocean 
migration  occurs  slightly  above  the  speed  at  minimum  t  (shown  with  an  arrow),  which  is  about 
\K,  .,(©)•  while  river  migration  takes  place  nearly  at  K,  .,(•)-  For  the  definition  of  Fcrlt  see  text 
and  the  caption  to  Fig.  9.  (Adapted  from  Brett'31, 3,)). 
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Theodore  Y.  Wn  and  George  T.  Yates 


Moss  (g) 

Fig.  9.  Metabolic  rate  of  sockeye  salmon  (Oncorhunchus  nerka )  as  a  function  of  body  mass  at 
various  levels  of  swimming  activity  (15'C).  The  standard  level  was  determined  by  extrapolation 
to  zero  swimming  speed  of  the  oxygen  consumption  vs.  speed  relation.  The  line  at  full  activity 
corresponds  to  the  maximum  speed  which  could  be  sustained  for  60  min,  Vtril.  Intermediate 
fraction  levels  correspond  to  swimming  speeds  equal  to  J,  and  J  of  n,.,  respectively. 


FIG.  3.  The  "metabolic  rate  for  swimming,"  based  on  meas¬ 
ured  oxygen  consumption  less  the  basal  rate,  of  sockeye  salmon 
fOncorhynchus  nerkaj  as  a  function  of  body  mass  at  various 
levels  of  activity  (15"C) ;  data  adapted  from  Brett  (1965a). 
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FIG.  Variation  of  swimming  speed  with  body  length  at 

specific  levels  of  activity .  Burst  speed  data:  ♦  dace ,  a 
goldfish ,  •  trout,  B+  barracuda ,  Ps+  porpoise ,  IH-  dolphin 

(Bainbridge,  1961).  Full  activity  data:  open  syml>ols 
( Bainbridge ,  1961);  ®  sockeye  salmon  (15°C)  (Brett,  1965a); 

0  bass,  ®  coho  salmon  (Dahl berg  et  al.#  1968);  % - - a  , 

C  cod,  R  red fish,  F  winter  flounder ,  S  sculpin,  P  pout  (Beamish, 
1966);  >-— <  goldfish  (Smit  et  al.,  1971);  e  larval  anchovy 
(tfunter,  1972);  » — «  herring  (Jones,  1963);  ®  salmon,  «  trout 
(Paulick  and  DeLacy,  1957).  1/4  -  max  data:  0  sockeye  salmon 

at  15°C  (Brett,  1965a).  The  lines  of  constant  Reynolds  number. 
Re,  and  lines  of  constant  specific  speed  (in  l  s~l)  are  indi¬ 
cated  for  reference. 
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FLUID-DYNAMIC  DRAG 


Kigtue  22.  Selected  results  on  the  drag  of  rotatio'nally- 
sMinneiric  bodies  (no  corrections  applied) . 
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FIG .  Variation  of  the  measured  metabolic  power  coeffic¬ 
ient,  Cp,  with  the  Reynolds  number  and  the  variation  of  the 
measured  dead-drag  coefficient ,  Cp^,  of  freshly  killed  salmon 
(adapted  from  Brett,  1963).  The  corresponding  values  of  the 
specific  speed  are  indicated  in  l  s~*  along  the  Cp  curve.  An 
estimate  of  the  swimming  drag  coefficient ,  Cp,  is  provided  by 
the  value  of  n Cp,  shown  with  a  typical  value  of  the  overall 
efficiency  n  -  0.25.  The  drag  coefficients  of  a  flat  plate 
for  both  laminar  and  turbulent  boundary  layers  are  shown  for 
reference. 
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FIG.  Scale  effects  in  the  variation  of  the  measured 

metabolic  power  coefficient  (for  swimming)  with  Reynolds 
number  for  five  different  size  groups  of  sockeye  salmon 
(adapted  from  Brett ,  1965a).  Solid  lines  are  least-square 
error  fit  to  Cp  at  specific  activity  levels.  Dashed  lines 
show  variation  of  Cp  with  various  activity  levels  for  each 
size  group .  Dash-dot  lines  illustrate  theoretical  thrust 
coefficient ,  CT,  computed  from  eg.  30  for  the  data  of 
Pyatetskiy  (1970b)  and  Webb  (1975). 
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FIG.  5.  Relation  between  specific  energy  cost  for  swlimlng, 
c,  and  body  mass  for  fish  swimming  at  different  levels  of 
activity.  Solid  lines  are  regression  lines  obtained  from  data 
of  Brett  (1965a)  for  salmon  at  1S°C. 


FIG.  6.  Comparison  of  the  observed  values  of  0  (Brett , 
1965a ,b)  and  y  (see  Fig.  5  and  Table  II)  with  the  similarity 
predictions  (eqs.  8,11)  based  on  the  three  distinct  reference 
states  characterized  by  a  laminar  boundary  layer,  a  turbulent 
boundary  layer,  and  CD  «  const,  (for  the  case  of  separated 
cross  flows)  . 
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Figure  13,  Dorsal  and  lateral  views  of  mackerel  (Scomber 
scomber)  showing  caudal  f inlets,  one  of  which  is  enlarged  to  show 
manner  in  which  a  layer  of  epithelium  attaches  the  f inlet  along 
most  of  its  length. 
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Figure  12.  Diagrammatic  section  through  integument  of  three 
fishes  which  have  subdermal  canals  and  associated  pores  directed 
obliquely  rearwards.  Ctenoid  scales  where  present  darkly  shaded. 
Direction  of  movement  of  the  fish  indicated  by  arrow,  a.  Ruvettus. 
b.  Tetragonous  (an  aberrant  stromateoid  where  the  scales  form 
underlying  channels),  c.  Schedolphilous.  Ctenoid  scales  absent. 
The  diagram  at  right  indicates  the  manner  in  which  the  convex 
leading  surface  of  the  fish  is  supposed  to  eject  sea  water  during 
the  swimming  cycle. 
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Fig.  11.  The  streamlines  (in  the  laboratory 
frame)  of  a  self-propelling  prolate  spheroid 
of  eccentricity  0.  9  (upper)  and  of  an  externally- 
driven  spheroid  of  the  same  shape  (lower 
figure,  with  v3  =  vn  =  0)- 


Fig.  12.  Streak  photographs  of  a  freely 
swimming  specimen  of  Paramecium 
enudatum  (upper)  and  of  a  dead  one  sedimen¬ 
ting  (towards  right)  under  gravity  (lower 
figure). 
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Hydromechanics  of  swimming  propulsion. 

Part  1.  Swimming  of  a  two-dimensional  flexible  plate  at 
variable  forward  speeds  in  an  inviscid  fluid 

By  T.  YAO-TSU  WU 

Californio.  Instituto  of  Technology, 

Pasadena,  Californio 

(Received  21  July  1970) 

The  most  effective  movements  of  swimming  aquatic  animals  of  almost  all  sizes 
appear  tohave  the  form  of  a  transverse  wave  progressing  along  the  body  from  head 
to  tail.  The  main  features  of  this  undulatory  mode  of  propulsion  arc  discussed  for 
the  case  of  largo  Reynolds  number,  based  on  the  principle  of  energy  conservation. 
The  general  problem  of  a  two-dimensional  flexible  plate,  swimming  at  arbitrary, 
unsteady  forward  speeds,  is  solved  by  applying  the  linearized  inviscid  flow  theory. 
The  largo-time  asymptotic  behaviour  of  an  initial-value  harmonic  motion  shows 
the  decay  of  the  transient  terms.  For  a  flexible  plate  starting  with  a  constant 
acceleration  from  at  rest,  the  small-time  solution  is  evaluated  and  tho  initial 
optimum  shape  is  determined  for  tho  maximum  thrust  under  conditions  of  fixed 
power  and  negligible  body  recoil. 


1.  Introduction 

Aquatic  animals  propelling  themselves  in  water,  or  in  other  liquid  media, 
span  a  wido  range  in  their  sizes  and  speeds.  Largo  cetaceans,  such  as  porpoises 
and  whales,  may  have  lengths  from  2  to  30  m,  and  can  swim  at  cruising  speeds  of 
from  6  to  10  m/s  (Lang  &  Pryor  1966).  Microscopic  organisms  such  as  paramecia 
and  spermatozoa,  ranging  from  300  (tm  down  to  50  fim  in  lciigth  with  length- 
diameter  ratio  from  20  to  100,  can  swim  at  speeds  from  1000  to  80/im/s.Inbctwccn 
these  two  extremities  there  are  many  species  of  fishes  and  aquatic  animals  of 
various  sizes.  Based  on  the  characteristic  length  l  of  a  body  moving  at  velocity 
V  in  a  liquid  of  kinematic  viscosity  v,  the  Reynolds  number,  It  =  Ul/v,  measures 
tho  relative  magnitude  of  tho  time  average  of  inertial  stress  to  viscous  stress. 
The  value' of  R  is  of  order  10s  for  the  most  rapid  cetaceans,  10“  for  migrating 
fishes,  105-10J  for  a  great  variety  of  fishes,  about  102  for  tadpoles,  down  to  about 
1  for  Turbalrix,  10-3  or  less  for  paramecia  and  spermatozoa  (Gray  1968,  p.  437), 
and  to  the  extreme  of  10-6  or  less  for  bacteria.  Thus,  the  Reynolds  number  R 
covers  practically  the  entiro  range  of  interest  known  to  hydi'odynamicists. 
Lighthill  (1969)  has  given  an  excellent  survey  of  tho  hydromechanics  of  aquatic 
animal  propulsion,  which  has  elucidated  both  the  zoological  and  hydromcchanical 
aspects  of  the  subject. 

Although  R  may  vary  greatly  from  case  to  case,  the  most  effective  movements 


of  swimming  propulsion  employod  by  a  largo  numbor  of  aquatio  animals  of 
drastically  different  sizos  havo  boon  obsorved  to  diflor  vory  littlo  from  an 
undulatory  motion  of  tho  body,  in  tho  form  of  a  transvorsc  wavo  propagating 
along  the  body  from  head  to  tail.  A  great  majority  of  many  species  of  fishes  can 
be  singled  out  as  a  pre-cminont  class  of  this  mode  of  propulsion.  Tho  rcmarkablo 
performance  of  some  cetaceans  (dolphin,  porpoises,  whales,  etc.)  and  some  well- 
known  gamo  fish  families  (tuna,  wahoo,  marlin,  swordfish,  etc.),  using  strong 
tails  of  large  aspoct  ratio,  is  only  a  variation  of  this  basic  undulatory  mode. 
In  tho  world  of  micro-organisms,  an  enormous  variety  of  croaturcs,  ranging 
from  minuto  bacteria,  larger  but  still  primitive  protozoa,  to  higher  level  sperma¬ 
tozoa,  havo  beon  observed  to  employ  cither  uniformly  propagating  transverso 
waves,  or  whip-liko  waves,  or  helical  waves  along  slender  flagella  as  principal 
means  of  propulsion.  The  basic  transverso  wavo  mode  thus  seems  to  bo  little 
affected  by  tho  Reynolds  number  over  such  a  wido  range.  Ilowovcr,  tho  funda¬ 
mental  principles  underlying  tho  hydromechanics  of  swimming  propulsion  do 
bccomo  very  different  for  large  or  small  values  of  the  Roynolds  number. 

For  Roynolds  number  large,  tho  swimming  propulsion  depends  primarily  on 
the  inertial  offect,  since  tho  flow  outside  a  thin  boundary  layer  noxt  to  the  body 
surface  is  irrotational.  Viscosity  of  tho  fluid  is  unimportant  except  in  its  role  of 
generating  tho  vorticity  shed  into  the  walco,  and  of  producing  a  thin  boundary 
layer,  and  hence  a  skin  friction  at  tho  body  surface.  As  the  body  performs  an 
undulatory  wavo  motion  and  attains  a  forward  momentum,  tho  propulsive 
force  pushes  tho  fluid  backward  with  a  net  total  momentum  equal  and  opposite 
to  that  of  the  action,  while  tho  frictional  rcsistanco  of  the  body  gives  riso  to  a 
forward  momentum  of  tho  fluid  by  entraining  soino  of  tho  fluid  surrounding  tho 
body.  Tho  momontum  of  roaction  to  tho  inertial  forces  is  concentrated  in  tho 
vortox  wako  duo  to  tho  small  thickness  and  amplitudo  of  tho  undulatory  trailing 
vortox  slicot;  this  backward  jet  of  fluid  expelled  from  tho  body  can,  however,  be 
counterbalanced  by  the  momentum  in  responso  to  tho  viscous  drag.  When  a  self- 
propelled  bodyis  cruisingataconstantspeed,thefonvardand  backward  momenta 
exactly  balance;  they  can  nevertheless  bo  evaluated  separately.  This  mechanism 
of  swimming  motion  at  largo  Reynolds  numbers  has  been  elucidated  by  von 
Karm&n  &  Burgers  (1943)  for  the  simple  case  of  a  rigid  plate  in  transverso  oscilla¬ 
tion.  Swimming  of  slender  fish  has  been  treated  by  Lighthill  (I960);  and  tho  waving 
motion  of  a  two-dimensional  floxiblo  plate  has  been  calculated  by  Wu  (1961). 

In  the  other  extremity,  movements  of  microscopic  bodies  always  correspond 
to  small  Reynolds  numbers.  The  propulsion  in  tills  range  depends  almost  ontirely 
on  tho  viscous  stresses,  since  the  inertial  forces  arc  then  oxtrcmcly  small,  except 
possibly  for  the  motions  at  very  high  frequencies.  Oscillatory  motions  in  a 
viscous  fluid  were  discussed  as  early  as  1801  by  Stokes.  Various  studies  of  tho 
swimming  of  microscopic  organisms  have  been  led  by  Taylor  (1951,  1952  a,  6), 
who  discussed  the  propulsion  of  a  propagating,  monochromatic,  transverso 
wavo  along  a  sheet  immersed  in  a  very  viscous  fluid,  and  later  evaluated  tho 
action  of  waving  cylindrical  tails  of  microscopic  organisms.  Further  studios  in 
this  field  havo  been  contributed  by  Hancock  (1953),  Gray  &  Hancock  (1955), 
Reynolds  (1965)  and  Tuck  (1968). 


1.3-34 


Fluid  Mech.  (1971),  vol.  '1ft,  part  3,  G21-S4-I 

Printed  in  Great  llritain 


521 


T.  Y.  I  Yu 


Hydromechanics  of  swimming  propulsion.  Part  2. 
Some  optimum  shape  problems 

lly  T.  YAO-T.SU  WU 

Cnliforniu  InRtiltito  of  Technology 
PoMmlonii,  Colifornio 

(Koroivotl  21  July  1970) 

The  optimum  shape  problems  considered  in  this  part  are  for  those  profiles  of  a 
two-dimensional  flexible  plate  in  time-hnrmonic  motion  that  will  minimize  the 
energy  loss  under  the  condition  of  fixed  thrust,  and  possibly  also  under  other 
isoperinietric  constraints.  First,  the  optimum  movement  of  a  rigid  plate  is 
completely  determined;  it  is  necessary  first  to  reduce  the  original  singular 
quadratic  form  representing  the  energy  loss  to  a  regular  one  of  a  lower  order, 
which  is  then  tractable  by  usual  variational  methods.  A  favourable  range  of  the 
reduced  frequency  is  found  in  which  the  thrust  contribution  coming  from  the 
leading-edgo  suction  is  as  small  as  possible  under  the  prescribed  conditions, 
outside  of  which  this  contribution  becomes  so  large,  as  to  be  bard  to  realize  in 
practice  without  stalling.  This  optimum  solution  is  compared  with  the  recent 
theory  of  Lighthill  (1970);  these  independently  arrived -at  conclusions  nre 
found  Lo  bo  virtually  in  agreement. 

Tlio  present  theory  is  further  applied  to  predict  the  movement  of  a  porpoise 
tail  of  large  aspect-ratio  and  is  found  in  satisfactory  agreement  with  the  experi¬ 
mental  measurements.  A  qualitative  discussion  of  the  wing  movement  in 
flapping  flight  of  birds  is  also  given  on  the  basis  of  optimum  eflicicnc}'. 

Tho  optimum  shape  of  a  flexible  plate  is  analysed  for  the  most  general  ease 
of  infinite  degrees  of  freedom.  It  is  shown  that  the  solution  can  be  determined  to 
a  certain  extent,  but  the  exact  shape  is  not  always  uniquely  determinate. 


1.  Introduction 

One  of  the  most  inspiring  questions  concerning  the  phenomena  of  aquatic 
animal  propulsion  and  of  flapping  flights  of  birds  and  insects  is  invariably 
connected  with  tho  highest  possible  hydrodynamic  efficiency.  This  problem 
has  been  brought  up  from  time  to  time  by  various  observers  who  have  noted  the 
impressive  capabilitity  of  these  animals  in  generating  fast  movements  at  low 
energy  cost.  According  to  the  first  principle  of  energy  balance  or  momentum 
consideration,  as  has  been  explained  in  part  1  of  this  paper  (Wu  1971),  much 
can  already  bo  said  about  the  desirable  shapes  of  body  movement:  that  nt  large 
Reynolds  numbers,  a  thin  two-dimensional  plato  gains  thrust  by  sending  a 
transverso  wave  from  head  to  tail,  with  amplitude  slightly  increasing  towards 
the  rear,  thereby  achieving  a  forward  swimming  velocity  somewhat  less  than 
the  phase  velocity  of  the  hod}'  wave  form.  As  for  tlu:  tail  of  large  aspect,  ratio 
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of  some  high-performance  fish,  the  tail  should  move  nearly  tangentially  lo  the 
path  traversed  in  the  space  by  the  body  wave  form.  These  basic  features  have 
been  elegantly  elucidated,  with  perhaps  more  physical  reasoning,  in  an  excellent 
roview  by  Lighthill  (1909).  However,  it  would  still  bo  of  great  interest  to  resol  vo 
a  quantitative  determination  of  the  optimum  shape  under  some  appropriate 
constraints. 

J  he  problem  of  tho  optimum  slmpe  is  interesting  in  its  own  right  from  tho 
mathematical  point  of  view,  sinco  the  effective  methods  of  solution  do  not 
seem  to  fall  into  the  known  categories  of  the  calculus  of  variation.  Tho  special 
case  of  a  two-dimensional  waving  plate  in  harmonic  motion  has  been  treated 
by  Wang  (1900),  who  adopted  a  discretized  Fourier  representation  of  the  body 
motion,  and  found  that  his  solutions  exist  only  for  a  set  of  eigenvalues.  How¬ 
ever,  it  is  found  in  tile  present  study  that  this  optimum  shape  problem  is 
bnsically  not  an  eigenvalue  problem,  and  therefore  merits  u  new  discussion. 
On  physical  grounds,  it  would  be  indeed  difficult;  to  see  the  significance  of  tho 
idea  that  the  shape  function  can  have  eigeusolutions. 

In  Ibis  part  wc  shall  consider  the  optimum  shape  problem  only  for  tho  case 
of  two-dimensional  flexible  plate,  of  negligible  thickness,  in  harmonic  motion. 
(Some  three-dimensional  problems  will  be  treated  in  part  3  of  this  paper.) 
Tho  two-dimensional  theory  is  reckoned  to  have  utility  in  problems  of  lifting 
surfaces  of  large  aspect-ratio,  such  ns  the  tails  of  some  cetaceans  and  high- 
performance  game  fish  (tho  lunate  tails:  swordfish,  tuna,  alhncore,  porpoises, 
etc.),  and  even  tho  wings  of  most  birds  and  some  insects.  The  optimum  shape 
problem  is  concerned  with  those  profiles  or  movements  that,  will  minimize  the 
energy  loss  under  the  condition  of  fixed  thrust  (required  to  overcome  the  viscous 
drag),  and  possibly  also  under  other  isoperinietric  constraints.  First,  the  optimum 
movement  of  a  rigid  plate  is  determined  b3r  reducing  the  original  singular 
problem  t.o  a  regular  one  of  a  lower  rank.  This  optimum  solution  is  found  to 
depend  on  two  variables:  one  being  the  reduced  frequency  and  the  other  a 
‘ proportional-loading  parameter’,  defined  ns  the  prescribed  thrust  coefficient 
divided  by  the  dimensionless  heaving  amplitude  squared.  For  given  loading 
parameter,  a  favourable  range  of  the  reduced  frequency  is  found  in  which  tho 
thrust  contribution  coming  from  tho  leading  edge  suction  is  ns  small  ns  possible 
under  the  prescribed  conditions.  This  consideration  seems  to  provide  the 
optimum  range  of  the  reduced  frequency  utilized  in  practice. 

These  theoretical  results  nre  further  applied  to  predict  the  movement  of  a 
porpoise  tail,  and  comparisons  made  with  the  experimental  investigation  of 
Lang  &  Daybell  (1909).  As  a  related  problem  of  interest,  the  optimum  movement 
of  a  flapping  wing  of  some  birds  or  flatfish  is  discussed  qualitatively. 

The  general  problem  of  optimum  shape  of  a  flexible  surface  having  an  infinite 
degree  of  freedom  is  finally  analysed  and  discussed.  It.  is  found  that  the  solution 
can  be  determined  to  a  certain  extent,  and,  with  the  additional  degrees  of  freedom, 
the  optimum  efficiency  can  be  further  improved  from  the  rigid-plato  value,  lmt 
tlie  exact  shape  is  not  uniquely  determinate. 
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Hydromechanics  of  swimming  propulsion. 

Part  3.  Swimming  and  optimum  movements  of 
slender  fish  with  side  fins 

By  T.  YAO-TSU  WU 

California  Institute  of  Technology, 

Pasadena,  California 

(Rocoivod  21  July  1070) 

This  paper  seeks  to  evaluate  the  swimming  flow  around  a  typical  slender  fish 
whose  transverse  cross-section  to  the  rear  of  its  maximum  span  section  is  of 
a  lenticular  shape  with  pointed  edges,  such  as  those  of  spiny  fins,  so  that  these 
side  edges  are  sharp  trailing  edges,  from  which  an  oscillating  vortex  sheet  is 
shed  to  trail  the  body  in  swimming.  The  additional  feature  of  shedding  of 
vortex  sheet  makes  this  problem  a  moderate  generalization  of  tho  paper  on  the 
swimming  of  slender  fish  treated  by  Lighthill  (19G0n).  It  is  found  here  that  tho 
propulsive  thrust  depends  not  only  on  the  virtual  mass  of  the  tail-end  section, 
but  also  on  an  integral  effect  of  variations  of  the  virtual  mass  along  tho  entire 
body  segment  containing  the  trailing  side  edges,  and  that  this  latter  effect  can 
greatly  enhance  the  thrust-making. 

The  optimum  shape  problem  considered  here  is  to  determine  the  transverse 
oscillatory  movements  a  slender  fish  can  make  which  will  produce  a  prescribed 
thrust,  so  as  to  overcome  the  frictional  drag,  at  the  expense  of  tho  minimum 
work  done  in  maintaining  the  motion.  The  solution  is  for  tho  fish  to  send  a 
wave  down  its  body  at  a  phase  velocity  c  somewhat  greater  than  the  desired 
swimming  speed  U ,  with  an  amplitude  nearly  uniform  from  the  maximum  span 
section  to  the  tail.  Both  the  ratio  I//c  and  the  optimum  efficiency  arc  found  to 
depend  upon  two  parameters:  the  reduced  wave  frequency  and  a  ‘proportional- 
loading  parameter’,  the  latter  being  proportional  to  the  thrust  coefficient  and 
to  the  inverse  square  of  the  wave  amplitude.  The  basic  mechanism  of  swimming 
is  examined  in  the  light  of  the  principle  of  action  and  reaction  by  studying  the 
vortex  wake  generated  by  the  optimum  movement. 


1.  Introduction 

Lighthill  (19G0n)  investigated  the  inviscid  flow  around  a  slender  fish  which 
makes  swimming  movements  in  a  direction  transverse  to  its  direction  of  loco¬ 
motion,  whilo  its  cross-section  varies  along  it  only  gradually.  Based  on  the 
slender-body  theory,  Lighthill  obtained  the  result  of  thrust  produced  by  the 
fish,  time-rate  of  work  clone  by  it,  and  the  rate  of  shedding  of  energy,  showing 
that  the  mean  values  of  these  quantities  all  depend  on  the  movement  and  body 
shape  at  the  tail-end  section  only,  and  that  they  will  vanish  with  the  virtual 
mass  of  the  tail.  What  has  primarily  been  implied  here  is  that  the  body  cross- 
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section  varies  so  gradually,  and  its  shape  is  so  smooth  (no  sharp  edges),  that 
the  cross-flow  remains  attached  to  body,  leaving  no  vortex  sheet  until  the  tail- 
end  section  is  reached.  This  situation  may  represent,  quite  accurately,  several 
wide  classes  of  aquatic  animals,  such  as  the  eel  Anguilla  vulgaris,  eel-like  fishes 
of  the  order  Heteromi,  ribbon  fishes  of  the  order  Allotriognathi,  members  of 
tho  order  Anacanthini  like  cod  and  similar  fishes,  as  have  been  discussed  in  an 
interesting  review  article  by  Lighthill  (1969). 

Thero  exist,  however,  other  classes  of  fishes  of  more  advanced  orders  which 
are  known  as  strong,  active  swimmers,  capable  of  putting  on  impressive  per¬ 
formances.  This  group  of  fish  orders  contains  members  of  the  order  Isospondyli 
(such  as  salmon,  trout,  etc.)  and  thoso  of  the  large  order  Ostariophysi,  to  which 
belong  most  of  the  successful  freshwater  fishes.  As  an  indicative  measure  of  the 
over-all  performance,  the  swimming  speed  U  is  usually  expressed  by  biologists 
in  units  of  body-length/scc,  or  n  =  U/l.  For  instance,  pike  (Esox  lucius),  salmon 
(Saltno  salar),  trout  ( Salmo  gairdneri),  dace  (Le.uciscus  leuciscus),  and  herring 
( Clupea  harengus)  have  been  observed  to  range  from  n  =  6/sec  to  12/sec  (see 
Hertel  1963,  chap.  G;  Gray  1968,  chap.  3).  Speeds  considerably  greater  than 
10/sec,  as  high  as  20/scc,  have  also  been  recorded  by  Walters  Sc  Ficrsteino  (19G4) 
for  tunny  ( Thunnus  albacores)  and  wahoo  in  open  water.  By  and  large,  the 
slenderness  parameter,  S,  defined  as  the  ratio  of  maximum  depth  (or  height,  or 
called  span  as  in  the  wing  theory)  to  body  length,  of  this  general  group  of  fishes 
is  moderate  to  small  (about  from  0-4  down  to  0-18  based  on  body  depth  including 
dorsal  and  ventral  fins  extended,  the  smallest  $  being  represented  by  the  fast- 
Bwimming  Pacific  saury,  Cololabis  saira ).  Their  thickness  in  the  third  dimension 
is  yet  smaller  than  their  span;  the  ratio  of  minor  to  major  axis  of  cross-section 
of  some  slightly  compressed  species  may  be  as  much  as  0-6  all  the  way  back  to 
the  caudal  pedunclo  (point  of  minimum  depth).  As  a  general  feature,  the  trans¬ 
verse  cross-sections  of  these  fishes  have  rather  rounded  edges  anterior  to  the 
section  of  maximum  depth,  turning  to  a  more  or  less  lenticular  shape  with 
fairly  pointed  edges  to  tho  rear  part  of  body  in  which  may  be  found  a  great 
variety  of  dorsal,  ventral,  pectoral,  anal,  and  possibly  other  smaller  fins,  to  be 
followed  by  tho  caudal  fin. 

As  for  the  detailed  fin  shapes  and  locations,  there  are  perhaps  as  many  different 
configurations  as  the  number  of  species.  However,  as  a  crudo  classification,  the 
dorsal  and  ventral  fins  may  be  arranged  for  hydromechanical  reasons  into  two 
main  types:  the  elongated  ‘ribbon  fin’  and  the  triangular  or  trapezoidal  ‘sail¬ 
shaped  fin’.  Some  pre-eminent  families  of  fast  fishes  equipped  with  conspicuous 
ribbon-fins  arc  dolphin  fish  Coryphaena  hippnrus,  which  has  a  dorsal  ribbon-fin 
all  the  way  from  head  to  caudal  peduncle  and  a  shorter  ventro-anal  one,  yellow- 
tail  Seriola  quinqueradiala,  atka  mackerel  Plcurogrammus  azonus,  porgies  and 
breams  (e.g.  Pagrus  major ,  Dcnlex  tumifrons).  Examples  of  fast  fishes  having 
sail-shaped  dorsal  and  ventral  fins  arc  mullet  Mugil  cephalus,  which  is  very  fast 
and  migrates  widely,  barracuda  Sphyraena  pi?iguis,  trout,  salmon,  nnd  herring. 
Also,  many  families  of  fishes  have  both  types  of  fins,  notably  bluefin  tuna 
Thunnus  thynnus,  skipjack,  and  some  mackerels.  Furthermore,  several  fishes 
(wahoo,  tuna,  skipjack,  mackerel  and  saury)  have  behind  the  main  fins  a  scries 
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OPTIMIZATION  PROBLEMS  IN  HYDROFOIL  PROPULSION* 

Th.  Yao-tsu  Wu,  Allen  T.  Chwang 
California  Institute  of  Technology,  Pasadena,  California 
and  Paul  K.  C.  Wang 
University  of  California  at  Los  Angeles 


This  paper  attempts  to  apply  the  principle  of  control  theory  to  investigate  the 
possibility  of  extracting  flow  energy  from  a  fluid  medium  by  a  flexible  hydrofoil 
moving  through  a  gravity  wave  in  water,  or  by  an  airfoil  in  gust.  The  present 
optimization  consideration  has  led  to  the  finding  that  although  the  flexible  hydrofoil 
may  have  an  infinite  number  of  degrees  of  freedom,  the  optimum  shape  problem  is 
nevertheless  a  finite -dimensional  one.  The  optimum  shape  sought  here  is  the  one 
which  minimizes  the  required  power  subject  to  the  constraint  of  fixed  thrust.  A 
primary  step  towards  the  solution  is  to  reduce  the  problem  to  one  of  minimizing  a 
finite  quadratic  form;  after  this  reduction  the  solution  is  determined  by  the  method 
of  variational  calculation  of  parameters.  It  is  found  that  energy  extraction  is 
impossible  if  the  incident  flow  is  uniform,  and  may  be  possible  when  the  primary 
flow  contains  a  wave  component  having  a  longitudinal  distribution  of  the  velocity 
component  normal  to  both  the  mean  direction  of  flight  and  the  wing  span.  When 
such  waves  of  sufficiently  large  amplitude  are  present,  not  only  flow  energy  but 
also  a  net  mechanical  power  can  be  extracted  from  the  surrounding  flow. 


*  This  paper  includes  further  extension  to  that  which  was  originally  presented  at 
the  Symposium, 

Lecture  Notes  in  Physics  21.  Optimization  and  Stability  Problems 
in  Continuum  Mechanics.  Edited  by  P.  K.  C.  Wang.  Springer- 
Verlag,  Berlin/New  York  (1973). 


I,  Introduction 

Some  previous  observations  on  fish  swimming  and  bird  flight  seem  to  suggest 
that  some  species  may  have  learned,  through  experience,  to  acquire  the  key  to 
high  performance  by  executing  the  optimum  movement  that  may  be  of  great  interest 
to  control  theory  related  to  fluid  mechanics.  An  especially  intriguing  aspect  of 
the  optimization  problem  concerns  with  the  possibility  of  extracting  energy  from 
surrounding  flow  by  an  oscillating  lifting  surface  (  such  as  the  fish  body  and  fins, 
bird  wings,  and  artificial  wings  like  airfoil  and  hydrofoil  )  and  its  associated 
effect  on  the  control  of  motion. 

This  general  problem  has  been  explored  to  various  degrees  of  generality. 

Based  on  the  approximation  of  potential  flow  with  small  amplitude,  it  has  been 
found  by  Ligh thill  (I960)  for  slender  bodies,  and  by  Wu  (1961  )  for  two-dimensional 
plates,  that  if  the  basic  flow  is  uniform,  energy  is  always  imparted  by  an  oscillat¬ 
ing  wing  to  the  surrounding  fluid,  and  an  extraneous  mechanical  work  must  there¬ 
fore  be  continuously  supplied  to  maintain  the  motion.  Even  though  it  is  impossible 
in  this  case  to  extract  energy  from  the  flow  field,  the  highest  possible  hydro¬ 
mechanical  efficiency  that  can  be  attained  by  a  wing,  subject  to  delivering  a  given 
forward  thrust,  can  be  very  high,  as  found  by  Wu  (1971  b,c)  for  the  two-dimen¬ 
sional  plate  and  a  slender  lifting  surface. 

As  was  subsequently  pointed  out  by  Wu  (1  972),  the  situation  becomes  drastical¬ 
ly  different  when  the  basic  flow  Is  no  longer  uniform,  but  contains  a  wave  compo- 
nent,such  as  gravity  waves  in  water,  or  wavy  gust  in  air.  The  contention  that  the 
wave  energy  stored  in  a  fluid  medium  can  be  utilized  to  assist  propulsion  has  been 
suggested  by  intuitive  observations.  Sea  gulls  and  pelicans  have  been  observed 
to  skim  ocean  waves  over  a  long  distance  without  making  noticeable  flapping 
motions  (save  some  gentle  twisting)  of  their  wings.  In  an  extensive  study  of  the 
migrating  salmon,  Osborne  (i  960)  found  that  the  increased  flow  rate  in  a  swollen 
river  did  not  slow  the  salmon  down  (for  known  biochemical  energy  expended  during 
the  travel)  by  that  much  a  margin  as  would  be  predicted  by  the  law  of  resistance 
in  proportion  to  the  square  of  their  velocity  relative  to  the  flowing  water.  Several 
possible  explanations  were  conjectured  by  Osborne,  including  the  prospect  that 
the  flow  energy  associated  with  the  eddies  in  river  could  be  converted  to  generate 
thrust.  To  explore  this  possibility  Wu  (1972)  introduced  an  energy  consideration 
to  an  earlier  study  of  Weinblum  (1954)  on  the  problem  of  heaving  and  pitching 
of  a  rigid  hydrofoil  In  regular  water  waves.  It  was  found  that  the  greatest 
possible  rate  of  energy  extraction  is  provided  by  the  optimum  mode  of  heaving 
and  pitching.  When  waves  of  sufficiently  large  amplitude  are  present,  not  only 
flow  energy  but  also  a  net  mechanical  power  can  be  extracted  from  the  wave 
field. 
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A  BIDIRECTIONAL  LONG- WAVE  MODEL 


Theodore  Yao-tsu  Wu 


Abstract.  For  modeling  weakly  nonlinear  and  weakly  dispersive  long  gravity 
waves  of  typical  amplitude  a  and  typical  length  A,  propagating  in  both  direc¬ 
tions  in  a  straight,  gradually  varying  channel  of  breath  b( x)  and  mean  water 
depth  h(x),  the  Boussinesq  equations  provide  a  versatile  model,  with  its  valid¬ 
ity  based  on  the  assumptions  that  a  =  a/h  <£  1,  e  =  (h/A)2  =  O(a)  and 
(d/dx)  logjb/i1/2)  =  0(t3/2).  On  this  same  basis,  a  new  bidirectional  long¬ 
wave  model  is  derived  to  evaluate  the  cross-sectional  mean  surface  elevation 
£( x,t )  =  (+  +  <-  +0,  where  £+  and  denote,  respectively,  right-going  and 
left-going  waves,  both  of  O(e).  and  (i  is  a  term  of  0(e2)  representing  the  interac¬ 
tion  between  £+  and  C--  The  evolution  equations  obtained  for  (+  and  (-  exhibit 
extensions  of  the  Korteweg-de  Vries  equation  to  comprise  the  additional  effects  of 
slow  variations  in  the  admittance,  ( bh */2),  of  varying  channels  on  evolving  waves. 
Main  features  of  this  model  include:  (i)  For  wave-and-channel-wall  interactions, 
this  model  accounts  for  reflection  and  transmission  of  long  waves  in  varying  chan¬ 
nels  while  maintaining  both  mass  and  energy  conserved  adiabaticallv.  (ii)  For 
wave-wave  interactions,  head-on  collisions  between  right-  and  left-going  solitary 
waves  are  shown  to  gain  a  total  phase  shift  which  is  an  algebraic  function  of  the 
amplitudes  of  the  colliding  waves,  (iii)  For  forced  generation  of  nonlinear  waves 
in  varying  channels,  this  model  admits  weakly  resonant  disturbances  with  both 
right-going  and  left-going  components. 


1.  Introduction 

The  general  subject  of  nonlinear  and  dispersive  waves  evolving  in  non-uniform  media 
has  been  of  strong  interest  and  active  development,  the  thrust  being  motivated  by 
attempts  for  the  ultimate  generality.  For  nonlinear  long  gravity  waves,  of  typical 
amplitude  a  and  length  A,  propagating  in  a  straight,  gradually  varying  channel  of 
breadth  b( i)  and  cross-sectional  mean-water-depth  h( x),  the  various  models  of  the 
Boussinesq  family  are  based  on  the  assumptions  that 

a  =  a//i  <  1,  e  =  (/i/A)2  =  0(a),  {d/dx)\og{b  h?)  =  0(e3/2),  (1) 

(see,  e.g.  Peregrine  [12];  Whitham  [22];  Miles  [9,  10,  11];  Wu  [19,  20,  21]).  The  models 
in  this  class  have  been  found  by  Teng  and  Wu  [17,  18]  to  be  in  broad  agreement  with 
existing  experiment. 

Along  a  different  approach,  the  additional  assumption  of  unidirectional  wave  motion 
in  varying  channels  has  been  applied  to  derive  the  Korteweg-de  Vries  (KdV)  class  of 
equations  by  Shuto  [15],  Miles  [11],  Teng  and  Wu  [17,  18]  and  others.  Similarly, 
variable-depth  form  of  the  Kadomtsev-Petviashvili  (K-P)  equation  has  been  given  by 
Liu  et  al.  [7].  Whilst  Shuto’s  equation  and  similar  models  are  known  to  be  adiabatic 
in  conserving  energy,  they  have,  however,  a  crucial  deficiency  in  not  conserving  mass, 
with  an  error  which  is  first  order  and  has  a  cumulative  effect.  This  effect  has  been 
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CONTROL  OF  THE  TURBULENT 
BOUNDARY  LAYER 

Professor  John  L.  Lumley 

CORNELL  UNIVERSITY 


The  construction  of  the  low-dimensional  models  of  Aubrey  et 
al.  is  briefly  outlined.  Motivated  by  a  desire  to  control  the  turbulent 
boundary  layer  via  feedback,  we  examine  the  effect  on  a  naturally 
occurring  streamwise  vortex  pair,  of  a  perturbing  pair  produced  by 
an  actuator.  We  show  how  the  resulting  velocity  field  modifies 
coefficients  in  the  above  low-dimensional  models  and  sketch  a 
control  strategy  which  reduces  the  bursting  rate  and  hence  the 
turbulent  drag  in  the  model  problem. 
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The  seminar  by  Professor  Lumley  was  canceled  at  the  last  minute  due  to 
the  sudden  illness  of  his  wife. 
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Situations  where  an  effectively  irrotational  freestream  contains 

regions  of  concentrated  vorticity  are  common  in  external  aerodynamics. 
Vortices  may  arise  as  a  consequence  of  shedding  from  some  upstream 

surface  (e.  g.  dynamic  stall)  or  near  certain  three-dimensional 
boundary  geometries  which  act  to  promote  vortex  formation.  Geometry- 
induced  creation  can  occur  in  any  situation  where  a  flow  along  a  wall 
approaches  a  surface-mounted  obstacle,  such  as  at  wing/hody  junctions 
and  near  computer  chips  mounted  on  electrical  circuit  boards.  Similar 
geometries  are  encountered  in  a  variety  of  internal  flows,  such  as 
branching  pipes,  and  in  turbine  and  compressor  passages.  The  cited 
examples  involve  macroscopic  large-scale  motions,  but  important  vortex 
motions  are  also  prevalent  at  much  smaller  scale.  Concentrated  small 
vortices  of  the  hairpin  type  are  known  to  be  an  important  influence  in 
a  boundary  layer  undergoing  transition  and  appear  to  be  the  dominant 

feature  in  the  production  of  turbulence  near  solid  surfaces. 

In  this  presentation,  recent  experimental  and  theoretical  work  on 

the  influence  of  vortex  motion  close  to  a  solid  surface  will  be  discussed, 

with  particular  reference  to  the  dynamics  of  turbulent  boundary  layers. 
When  the  flow  is  at  high  Reynolds  numbers,  boundary  layers  near  a 
surface  are  found  to  respond  in  a  generic  way  to  the  motion  of  the 
vortices  above;  in  this  process  an  abrupt  and  sharply  focused  eruption 
of  the  near-wall  fluid  invariably  develops.  Such  an  event  may  be 
characterized  as  an  unsteady  separation  of  the  surface  layer;  it  leads  to 

a  strong  viscous-in  viscid  interaction  with  the  external  flow  and 
culminates  in  the  violent  ejection  of  concentrated  near-wall  vorticity 
away  from  the  surface.  Generally,  it  is  not  possible  to  calculate  the 
evolution  of  these  ejections  at  high  Reynolds  number  using 

conventional  numerical  approaches.  Recent  progress  utilizing 

Lagrangian  methods  will  be  discussed. 
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SOME  APPLICATIONS 


1 .  VORTICES  CREATED  FROM  SHREDDING  FROM  AN  UPSTREAM  SURFACE 

•  AIRCRAFT  TRAILING  VORTICES 

•  DYNAMIC  STALL  ON  AIRFOILS  AND  HELICOPTER  BLADES 

•  TURBINES 

2.  GEOMETRY-INDUCED  VORTEX  CREATION 

•  WING/BODY  INTERSECTIONS 

•  END-WALL  BOUNDARY  LAYERS  IN  TURBINES 

•  COMPONENTS  ON  CIRCUIT  BOARDS 


3 .  TURBULENCE  AND  TRANSITION 


TURBULENT  BOUNDARY  LAYERS 


OBJECTIVES 

1 .  DETERMINE  CAUSE  AND  EFFECT  RELATIONSHIPS 

•  RELATION  OF  FLOW  STRUCTURE  TO  OBSERVED  EVENTS 

(PRIMARILY  NEAR  THE  SURFACE) 

•  ORIGIN  OF  THE  "LOW-SPEED"  STREAKS 

•  CAUSE  OF  WALL-LAYER  BURSTING 

(LOCAL  WALL  LAYER  BREAKDOWN,  INTERACTION) 

2.  ESTABLISH  DYNAMICAL  PICTURE 

•  CONTROL  ISSUES? 

APPROACH 

•  UNDERSTAND  RELEVANT  BASIC  PHYSICAL  PROCESSES 

•  IDENTIFY  MOST  BASIC  ELEMENT  OF  TURBULENCE  NEAR  A  SURFACE 

•  ANALYSIS  IS  ASYMPTOTIC  FOR  LARGE  RE 
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(a)  Symmetric 


(b)  Asymmetric 


HAIRPIN  VORTEX 

UNDERSTAND  INTERACTIONS  WITH 

(1)  THE  SURFACE  FLOW 

(2)  THE  BACKGROUND  SHEAR  FLOW 

(3)  EACH  OTHER 

THEN  IT  IS  POSSIBLE  TO  INTERPRET  OBSERVED  BEHAVIOR  IN  TRANSITION 
AND  THE  TURBULENT  WALL  LAYER 
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(a)  Symmetric 


(b)  Asymmetric 


HAIRPIN  VORTEX 

UNDERSTAND  INTERACTIONS  WITH 

(1)  THE  SURFACE  FLOW 

(2)  THE  BACKGROUND  SHEAR  FLOW 

(3)  EACH  OTHER 

THEN  IT  IS  POSSIBLE  TO  INTERPRET  OBSERVED  BEHAVIOR  IN  TRANSITION 
AND  THE  TURBULENT  WALL  LAYER 
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REGENERATION 


•  THE  PRODUCTION  OF  NEW  TURBULENCE 

•  RELEVANT  EVENTS  IN  TRANSITION  AND  TURBULENT  BOUNDARY 
LAYERS  ARE  UNSTEADY  SEPARATION  OF  THE  SURFACE  LAYER 

VORTEX  MOTION  PROVOKES  SURFACE  LAYER  ERUPTIONS  =>  NEW 
HAIRPIN  VORTICES 


BENEFITS 


1 .  MEAN-FLOW  PREDICTION  (WALL  LAYER) 

ATTACHED  FLOW  ON  EXTENDED  SURFACES 

=>  ALGORITHMS  FOR  (A)  HEAT  TRANSFER 

(B)  SUPERSONIC  BOUNDARY  LAYERS 

(C)  3-D  BOUNDARY  LAYERS 


2.  CONTROL 

•  ENHANCEMENT  (HEAT  TRANSFER,  MIXING) 

•  SUPPRESSION  (DRAG  AND  NOISE  REDUCTION) 
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SUMMARY 


1 .  GENERIC  ERUPTIVE  PHENOMENON 
(VORTEX-INDUCED  SURFACE-LAYER  ERUPTIONS) 

=>  VORTICITY  CONCENTRATIONS  NEAR  SURFACE 

=>  AN  ABRUPT,  DISCRETE  ERUPTION 

2.  VORTICITY  FIELD  ABOVE  THE  SURFACE 

=>  GROWTH  TO  LARGER  SCALE 

3.  WALL-LAYER  STREAKS 

4.  REGENERATION  OF  HAIRPIN  VORTICES  NEAR  THE  SURFACE 
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VORTEX  INTERACTIONS  WITH  WALLS 


•  ALL  VORTICES  NEAR  A  SURFACE  INDUCE  AN  ADVERSE 
PRESSURE  GRADIENT  IN  THE  NEAR-WALL  FLOW 


•  AT  HIGH  RE,  THIS  PROVOKES  SEPARATION  OF  THE  SURFACE 
LAYER 


=>  A  GENERIC  EVENT  IN  TWO  (OR  THREE)  DIMENSIONS 


VORTICITY  NEAR  THE  SURFACE  IS  ABRUPTLY  COMPRESSED 
ALONG  A  NARROW  STREAMWISE  BAND  AND  LEAVES  THE  SURFACE 
IN  A  SHARPLY-FOCUSED  SPIKE  (DOUBLE-SIDED  SHEAR  LAYER) 


A  VISCOUS-INVISCID  INTERACTION  IN  WHICH  "PARCELS"  OF 
VORTICITY  ARE  EJECTED  INTO  THE  EXTERNAL  FLOW 


General  Sequence  of  Events 


Stage  1  -  Triggering  Phase 

An  adverse  pressure  gradient  leads  to  the  evolution  of 
recirculation  in  surface  flow  (1st  MRS  condition) 

=►  zero  vorticity  line  (2-D) 
zero  vorticity  surface  (3-D) 


Stage  2  -  Separation 

Focussing  and  concentration  of  surface  layer  vorticity 
=*-  outward  moving  spike 


i 

# 


/* 


Onset  of  Interaction 

•  sharp  local 

pressure  changes 


developing  narrow-band 
eruption  (localized) 
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Stage  3  -  Strong  Interaction 


Double-sided  vorticity  layer  leaves  the  surface;  a  window 
opens  in  boundary  layer 


Stage  4  -  Inviscid  Interaction 


Last  three  stages  involve  moving 
fronts  of  concentrated  vorticity 
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cross-vortex  plane 


2(c) 


Figure  2.  Instantaneous  flow  induced  by  a  vortex  near  a  wall  in  a  frame  of  reference 
moving  with  the  vortex  core  and  in  the  cross- vortex  plane:  (a)  instantaneous 
streamline  patterns  (...  possible  spiral  motion),  (b)  flow  speed,  pressure,  and 
pressure  gradient  induced  near  the  surface,  (c)  the  cross-vortex  plane. 
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(a)  The  evolution  of  recirculation  (b)  Eddy  growth 

near  the  surface 


(c)  Focussing  and  compression;  solid  lines  denote  the 
equi-vorticity  contours  in  the  surface  layer;  arrows 
indicate  compression  of  the  erupting  spike 


Figure  3.  Schematic  diagram  of  the  stages  in  the  separation  of  a 
two-dimensional  surface  layer  (not  to  scale). 
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Figure  4.  Temporal  development  of  the  displacement  thickness;  plotted  curves 
are  at  t  =  0.25  (0. 10)  0.95  and  t.  =  0.989. 


Figure  5.  Schematic  of  boundaiy-Iayer  structure  near  a  point  of  eruption,  (a)  Structure  near 
x,  (not  to  scale);  (6)  a  typical  velocity  profile  near  x  =  x,  (for  a  stationary  wall). 


A  temporal  sequence  of  end-view  photographs  illustrating 
streamwise  necklace  vortex  legs  and  associated  surface  eruptions 


Figure  5.  End-view  of  necklace  vortex  legs  and  induced  surface  eruptions. 

Visualization  using  hydrogen  bubbles  and  light-sheet  illumination. 
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hairpin  streamwise 

vortex  vortex 


a)  End-view,  light-sheet  photograph  showing  the  breakdown  of  an  eruptive 
spire  (low-speed  ridge)  induced  by  a  streamwise  necklace  vortex  leg. 


process  for  the  low-speed,  eruptive  ridge. 


Figure  20.  Characteristics  of  the  breakdown  of  the  eruptive  spires  induced  by 
streamwise  necklace  vortex  legs. 
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(a) 


Head 


Subsidiary 


(c) 


Figure  6.  Evolution  of  a  symmetric  hairpin  vortex  in  a  shear  flow. 

(a)  initial  distortion,  (b)  development  of  vortex  legs  and  head, 
(c)  evolution  of  subsidiary  vortices  and  penetration  toward  the 
surface. 
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II 


2.1-28 


Figure  10.  Evolution  of  an  asymmetric  hairpin  vortex  in  a  turbulent  mean  profile:  (a) 
top  view,  (b)  side  view,  (c)  end  view;  c  =  0.0013,  400  time  steps  with 
At  =  0.016.  The  vortex  position  is  plotted  every  40  time  steps. 


2.1 


Figure  12.  Evolution,  of  multiple  hairpin  vortices  in  a  uniform  shear;  (a)  top  view,  (b) 
side  view,  (c)  end  view;  all  three  coordinates  are  plotted  on  the  same  scale. 
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SUBCRrTICAL 
LAMINAR  FLOW 


CONSTANT  LEVEL 


PIPETTE 

(SYRINGE) 


Figure  14.  Experimental  system  for  generating  single  hairpin  vortices  using  controlled  injection 
through  a  narrow  streamwise  slot  into  a  subcritical  laminar  boundary  layer. 


LEGS 


PLAN-VIEW 


FLOW  ►► 


SIDE-VIEW 


(a) 


Head 


Figure  15.  Illustration  of  experimental  generation  of  a  single  hairpin  vortex  by 

surface  injection  using  system  shown  in  figure  14.  (a)  Dual-view  picture 
of  dye-marked  single  hairpin  vortex,  (b)  Isometric  schematic  of  single 
hairpin  after  generation. 
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B 


BW 


Figure  17.  Dual-level,  plan-view  visualization  of  hydrogen  bubble  patterns  generated  by  left-to-right 

passage  of  developing  hairpin  vortex  (c.f.  figure  16).  BW  denotes  the  position  of  the  bubble  wire, 
P  the  location  of  an  initial  "pocket"  formation  (Falco,  1991),  L  the  location  of  the  vortex  legs,  S 
the  evolution  of  secondary  hairpin  vortices  (see  §5),  and  Cl  the  plane  of  symmetry.  Bubbles  are 

generated  at  y/8=0.4  for  upper  row  visualizations,  and  y/5=0.1  for  lower  row  visualizations. 


Plan-view  photograph  showing  a  streamwise  necklace  vortex 
leg  (visualized  by  upstream  hydrogen  bubble  wire)  interacting 
with  surface  fluid  (visualized  by  bubble  wire  near  surface)  to 
create  low-speed  streak  at  surface. 


Figure  1 8  Plan-view  of  necklace  vortex  leg  interaction  to  create  low-speed  streak 


Primary  Hairpin 


Interaction 
with  Legs 


Trailing  Low-Speed 
Streak  Regions 


(a)  Symmetric  case  with 
wide  separation  of  legs 


Interaction 
with  Legs 


Trailing  Low-Speed 
Streak  Region 


(b)  Symmetric  case  with  legs 
closer  together;  the  two 
streaks  merge 


Single-Leg  ,, 

Hairpin  ~  ** 


Interaction 
with  Leg 


Trailing  Low-Speed 
Streak  Region 


(c)  Asymmetric  case 


Figure  19.  Schematic  diagram  of  the  processes  whereby  moving  hairpin  vortices 
induce  low-speed  streaks. 
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Hairpin  Leg 

Surface 
Layer 


Focussed  Eruptive 
Growth 


(a)  Symmetric  hairpin  legs 


Figure  21. 


Schematic  diagram  of  location  of  vortex-induced  separation  of  the 
surface  layer  for  symmetric  and  asymmetric  hairpin  vortex  legs. 
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b)  Rapid  Outward  Movement  of  the  Erupting 
Ridge  Which  Contains  Concentrated  Vorticity. 


e)  Complete  Generation  of 
Secondary  Hairpin  Vortices 


2.1-42 


FLOW  - 


Figure  23. 


HBW 


a)  t=0.0  Sec. 


b)  t=0.166 


c)  t=0.25 


B 

C  ! 


t=0.416 


Plan-view  hydrogen  bubble  wire  visualization  sequence  illustrating  the  development  of 
secondary  vortices  near  the  surface  as  a  primary  hairpin  vortex  passes  a  fixed  streamwise 
location.  HBW  denotes  the  position  of  the  hydrogen  bubble  wire,  HPH  is  the  location  of  head 
of  primary  vortex,  A  the  location  of  the  trailing  legs  of  primary  vortex,  B  the  development  of 
a  secondary  vortex  behind  head  of  primary,  C  the  development  of  secondary  vortices 
adjacent  to  the  legs  of  the  primary  vortex,  and  L1,L2  the  legs  nearest  the  symmetry  plane  for 
the  secondary  vortices  indicated  by  C. 
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Figure  24.  The  generation  of  secondary  vortex  via  surface  interaction  for  an  asymmetric  hairpin 


vortex 
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Hydrogen 


Light  Sheet 
Illumination  of  Bubble 
Cross-Section 


Hydrogen 

Bubble 

Sheet 


Laser 

Light 

Sheet 


Figure  25. 


Schematic  diagram  of  laser  sheet  visualization  technique  employed 
to  visualize  eruptive  events  in  a  turbulent  boundary  layer. 
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(a)  tip  penetration,  y+=92 


(b)  tip  penetration,  y+=105 


Figure  26.  End-view  visualization,  using  a  horizontal  hydrogen  bubble  wire  and  laser  light-sheet 

illumination  (see  figure  25),  showing  eruptive  spires  emanating  from  near-wall  region  of  a 
turbulent  boundary  layer  (Ree=1150).  Lateral  field  of  view  *  250  wall  units,  wire  located  y+  =7 
from  surface,  and  x+  =»  150  upstream  of  light  sheet.  Compare  with  eruptions  induced  by 
streamwise  vortex  as  shown  in  figure  5. 
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SYNTHETIC  TURBULENCE 

Professor  K.  R.  Sreenivasan 
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We  discuss  schemes  for  generating  stochastic  signals  that 
possess  many  properties  of  high-Reynolds-number  turbulence.  The 
use  of  such  signals,  designated  "synthetic  turbulence",  in  providing 
initial  conditions  for  direct  numerical  simulations  and  in  the 
development  of  turbulence  models  will  be  discussed  briefly. 
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SYNTHETIC  TURBULENCE 


References 

"Synthetic  turbulence"  by  A.  Juneja,  D.P.  Lathrop,  K.R . 
Sreenivasan  and  G.  Stolovitzky,  Phys.  Rev.  E.  May  1994 

"A  scheme  for  generating  initial  fields  for  the  direct 
numerical  simulation  of  box-turbulence",  by  A.  Juneja,  G . 
Erlebacher,  and  K.R.  Sreenivasan,  preprint,  April  1994 

Synthetic  turbulence  (ST)  =  turbulence-like  stochastic 
signals  constructed  from  some  simple  rational  "rules" 

Among  other  things:  ST  has  the  same  energy  spectrum,  dissipation 
rate,  single-point  PDFs  for  velocities  and  velocity  increments, 
intermittency  characteristics,  structure  functions  up  to  some  high 
order,  fractal  dimension,  etc.  as  high-Reynolds-number  turbulence 

Two  practical  goals 

Construct  ST  to  provide  efficient  intial  conditions  for  DNS 

Obtain  subgrid-scale  viscosity  and  a  few  other  modeling 
constants  (e.g.,  Kolmogorov  constant,  turbulent  Prandtl 
number,  some  constants  in  the  K-8  model) 

IF  THIS  IS  POSSIBLE,  THE  "RULES"  MAY  BE  THOUGHT  i 
TO  CONTAIN  THE  RIGHT  PHYSICS  i 


FOCUS  ON  PRINCIPLES  BY  CONCENTRATING  ON  ONE¬ 
DIMENSIONAL  TRACE  OF  VELOCITY 
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Generally,  initial  conditions  are  chosen  to: 

(a)  Satisfy  continuity  equation 

(b)  Possess  the  desired  energy  spectrum 

(in  particular,  kinetic  energy  and  energy  dissipation  rate) 

Shortcomings 

a.  For  homogeneous  and  isotropic  turbulence,  the  decay 
history  is  set  largely  by  the  initial  state 

b.  Slow  convergence  of  high-order  statistics  and  large-scale 
properties 

For  1283  box,  random  initial  conditions  take  about  2500 

iterations  for  convergence 
(each  step  taking  about  3  sec  on  Cray-YMP). 

cf:  about  300  iterations  if  started  with  ST 

(comparison) 

c.  Total  vorticity  increases  initially  before  decaying 
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Non-universal  and  quasi-universal  aspects 

(a  schematic) 


successive  break-up  of  eddies 


NOTATION 

8:  energy  dissipation  rate  per  unit  mass 
er:  average  of  e  over  a  box  of  size  r 

Aur  =  Au(r)  =  u(x+r)  -  u(x) 
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£  Ax 


Multiplicative  Process 


Basis  for  the  multifractal  behavior  discussed  in:  C.  Meneveau  & 
K.R.  Sreenivasan,  J.  Fluid  Mech.  224, 429-484  (1991) 
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^istK  buttons 


A  fluftSi-  McJUL 

choose 

P(M)  =  0.55[M-Mo]  +  0.58[M-(l-Mo)] 

and  determine  M0  by  requiring  that 

<Mq>m(Kjei  =  <M^>  exp  (1) 

for  some  values  of  q. 

q  =  0:  normalization  condition 

q  =  1:  distribution  average  =  0.5 

q  =  2:  <M2>model  =  0.5[Mo2+  (1-M02)] 

2 

Computing  ^eXp  from  previous  figure, 

we  get,  using  (1),  the  result  that 

M0  =  0.7 

"beKvokon.  of  tW.e  (a- model  from,  multiplier  distri button 
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Figure  3:  A  comparison  between  moments  computed  from  the  measured 
multiplier  distributions  and  those  computed  for  the  different  models 
considered  in  the  text.  Experimental  data  were  obtained  from  a  record 
length  of  810,000  data  points.  The  convergence  of  moments  was 
reasonable;  for  example,  in  the  last  half  decade  of  the  record  length,  the 
variations  observed  were  smaller  (in  the  log  scale)  than  the  symbol  size. 
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d(x) 


Step  2:  Obtain  Aur  from  er  using  KRSH 

A.N.  Kolmogorov,  J.  Fluid  Mech.,  12,  82 


Aur  =  V  (rer)l/3 

already  discussed 


V  is  a  stochastic  variable  independent  of  r  and  rer  and  is 
hence  ’universal’ 


1.  Is  V  indeed  universal? 

Appears  to  be  so. 

G.  Stolovitzky,  P.  Kailasnath  and  K.R.  Sreenivasan,  Phys.  Rev.  Lett. 
69,1178(1992) 

A.  Praskovsky,  Phys.  Fluids  A,  4,  2589  (1992) 

S.  Chen,  G.  Doolen,  R.  Kraichnan  and  Z.-S.  She,  Phys.  Fluids  A,  5, 
458  (1993) 

I.  Hosokawa,  J.  Phys.  Soc.  Japan,  62,  10  (1993) 

S.  Grossmann  and  D.  Lohse,  Phys.  Fluids  A,  6,  611  (1994) 

2.  Why  is  V  universal,  and  what  are  its  properties? 

G.  Stolovitzky  and  K.R.  Sreenivasan,  Rev.  Mod.  Phys.  66, 229-240  (1994) 

True  for  all  processes  resembling  "anti-persistent"  fractional  Brownian 
motion;  "conditioning"  weakens  the  connection  to  the  specificity  of  the 

process 
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log10<  I  Au|  |  (rer)1/3>  log10P(V|  (r£r)v 


FIG.  2, 


G.  Stolovitzky  and  K.R.  Sreenivasan,  Phys.  Rev. 
Lett.  69, 1178, 1992 

Atmospheric  data,  6  m  above  ground,  R^  between 
1800  and  2500 
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Step  3 

The  construction  of  the  velocity  signal 


I  step 


't 


aw,  -  v,  (re,), 


1^+  AUf 


II  step 


^  t  M, .  V,  Lrt,)* 
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(fcl  CEMKftAL) 

mismatch  at  1/2 


Uj/2  =  a(u0+Au|)  +  (l-a)(uj-Aii2) 

a  =  1/2:  random  additive  process 

(  by 

For  all  a,  can  show  that:  mrtK*^) 


,  fO  if  p  =  2m+l 
<(A'ur)^>-v  j  r 

(,  if  p  =  2m 


Need  "structure"  somehow 


Tertiary  cascades  provide  such  a  structure 
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Figure  G:  Basic  lent  functions  and  used  in 
models. 
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Figure  7:  Example  velocity  sections  from  (a)  atmospheric  turbulence,  (b) 
synthetic  turbulence  based  on  the  p-model,  (c)  synthetic  turbulence  based 
on  the  /-model. 


How  did  the  large  scale  behavior  appear? 


Simply  by  allowing  several  largest  inertial-range 

scales  to  co-exist 


Why  does  this  simple  scheme  MworkM? 

Because  of  the  existing  balance  between  the  large 
and  small  scales  set  by  the  right  dissipation  rate 


Would  it  work  in  other  flows  as  well? 

M  - ....  .1  . . .  I  I  »I|».W  ■  »■  ■■  ~ 

Probably  not  -  certainly  not  where  the  large  and 
small  scales  are  not  in  equilibrium 
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Synthetic  Turbulence  Model  in  3-Dimensions 

Generating  a  divergence-free  artificial  velocity  field 

Generate  3  scalar  fields  u',  v'  and  w' 

Define  uF  =  (uf,v',w'). 

uF  can  be  decomposed  as 
uF(x)  =  u(x)  +  uc(x) 
where 

V  •  u  —  0  ,  Vxuc  —  0. 


This  decomposition  can  be  accomplished 
in  Fourier  space  according  to 


k  uF 

k 2 


Transform  u  back  to  physical  space  to  define 
a  divergence-free  artificial  velocity  field. 


Effect  of  Projection  ? 
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The  Kolmogorov  constant  and  the 
intermittency  exponent 


4>UG0  =  Ck<e>2/3k~5/3 (kL)  jj 

Kolmogorov  constant  Intermittency  correction 


Au(r)  =  u(x+r)  -  u(x) 

known  from  analysis  and  construction 


can  therefore  obtain  the  quantity  (|)u(k),  thus  and  ft 


Ck  =  (2/tc)  r(l+C2)  sinK2/2) 
=0.51 

(cf:  best  experimental  estimate  ~  0.5) 


ii  =  C2  -2/3 
=  0.027 


(cf:  best  experimental  estimate  0.025±0.005) 


For  the  p-model, 
C2=log2{[0.72/3+0.32/3]/2} 
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Lftrge.-ed.dij  simulations  Omd  Sftkgrtd  see  t  Viscosity 

1.  Typical  scale  ratio 

3/4 

L/r|  =  Re^4 

2.  Subgrid  scale  models  for 
large-eddy  simulations 

a.  Resolve  only  scales  r  and  above 

b.  Incorporate  effects  of  unresolved 
scales  as  an  effective  viscosity 

c.  Rewrite  NS=0  in  terms  of  ur,  and  the 

2 

additional  terms  as  <ve(r)>V  ur 

3.  Computation  of  <ve(r)> 

a.  Write  <ve(r)>  =  a  <IAurl>  r 

b.  Use  Aur  =  V(rer)1/3 

c.  Use  p-model  or  alternatives 

d.  Impose  the  requirement  that 
<ve(r)>  -*  v  as  r  -*  viscous  scale 
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Some  results 


1 .  Subgrid  scale  viscosity 

,  N  1/3  4/3  ,  n  on-i/3 

<ve(r)>  =  0.039  <e>  r  (r/L)  1 
< t ve (r) ] q>/<v e (r) >  =  <IVIq>/<IVI>q  Req  (r/L)^+q, 


where  the  exponents  |i  are  known  exactly: 

u  =  -logb<[iq/3> 


2.  The  constant  C..  in  the  K-e  model 
As  r  approaches  the  large  scale,  v  (r) 
where  vT  =  0.084  K2/<e> 


v. 


3.  Others 

a.  Corrsin-Obukhov  constant,  Cq  =  2.2/C^ 
(turbulent  Prandtl  number  =  0.9) 

b.  Intermittency  correction  for  the  scalar,  |Tq  =  0.039 


c.  For  uniformly  sheared  flows 
Cn=0.084/(l+bS2), 
where  b  can  be  computed. 

d.  Log-law  velocity  profile  can  be  generated  for  channel 

flows 
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Concluding  remarks 

1.  Analysis  of  high-Reynolds-number  turbulence 
data  reveals  much  statistical  order,, 

2.  This  order  can  be  codified  into  some  very  simple 
"rules". 

3.  The  "rules"  can  in  turn  be  implemented  to 
generate  synthetic  turbulence,  which  resembles  real 
turbulence  to  a  remarkable  degree. 

(Suggests  that  the  abstraction  of  the  "rules"and  their 
implementation  for  creating  ST  are  substantially 
correct  at  the  level  considered.) 

4.  ST  can  be  used  for  prescribing  efficient  initial 
conditions  for  DNS,  leading  to  much  saving  in 
convergence  time. 

5.  ST  can  also  be  used  to  evaluate  some  turbulence 
modeling  constants. 

6.  Extension  to  scalars  is  quite  straightforward. 
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Appendix 


Biographical  Information  on  Seminar  Speakers 


CURRICULUM  VITAE 


SERGEI  KASHIN 


[Pll  Redacted] 

Research  Experience: 

1993  to  present  -  Adjunct  Professor  of  Biology,  Department  of 
Biology  and  Marine  Science  Center,  Northeastern  University  and  East  Point, 
Nahant,  Massachusetts.  Investigating  the  motor  behavior  and  physiology  of 
motor  systems  in  lower  animals.  Teaching  graduate  course  “Biology  of 
Fishes”. 

1991-1992  -  Research  Fellow,  Environmental  Group,  ARCO  Oil  and 
Gas  Company.  Investigated  toxicity  of  industrial  effluents  on  aquatic 
ecosystems. 

1989-1991  -  Visiting  Researcher.  Scripps  Institution  of 

Oceanography.  In  vivo  study  of  fish  sarcomeres. 

1971-1989  -  Senior  Researcher,  Institute  of  Oceanology,  Moscow, 

USSR.  Short  term  changes  in  ocean  ecosystems. 

1979-1989  -  Director,  Laboratory  of  Ecology  and  Modelling, 

institute  of  Oceanology,  Academy  of  Sciences  of  USSR. 

1985  -  Chief  Scientist.  Amazon  River  Expedition.  Fish  adaptations  to 
the  low  oxygen  conditions. 

1988  -  Chief  Scientist,  Expedition  to  Benguella  Upwelling  Area, 
Southeast  Atlantic.  Productivity  in  upwelling  ecosystems. 

EDUCATION:  BA  (Diploma  of  Zoology),  Department  of  Biology,  Moscow 
State  University,  Moscow,  USSR,  1962. 

Dissertation:  Biology  of  Siberian  Sturgeon,  Acipenser  baeri. 

Ph.D.  (Candidate  of  Biological  Sciences).  Institute  of  Oceanology, 
Academy  of  Sciences  of  USSR,  Moscow. 

Title:  Organization  of  Fish  Locomotion  and  Geometry  of  Musculature. 
RESEARCH  GRANTS  AND  FELLOWSHIPS: 

Karolinska  Institute  of  Neurophysiology,  Stockholm,  Sweden, 
Postdoctoral  Fellow.  1975-1976, 

British-Soviet  Exchange  of  Scientists,  Bristol,  United  Kingdom. 

1977. 

U.S.  Geological  Survey/University  of  California  at  Los  Angeles, 
Department  of  Biology,  Visiting  Scientist.  1979-1980. 
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CURRICULUM  VITAE 


NEWMAN,  Barry  George 

CURRENT  APPOINTMENT:  Canadair  Professor  of  Aerodynamics,  McGill  University 


EDUCATION: 


St.  John’s  College,  Cambridge,  Strathcona  Open  Major  Scholarship  1944 
in  Mathematics 

Part  1  Mechanical  Sciences  Tripos,  Cambridge;  1st  Class  1946 

Pari  II  (Aeronautics);  1st  Class;  B.A.  and  M.A.  1947 

Ph.D.,  Sydney  University,  Australia  1951 


PROFESSIONAL  APPOINTMENTS: 


Scientific  Officer,  Flight  Research  Projects  1951-1953 

Assistant  Research  Officer:  Flight  Research  1953-1955 

University  Lecturer  1955-1958 

Visiting  Lecturer  1957 

Canadair  Visiting  Professor  in  Fluid  Mechanics  1958 

Canadair  Professor  of  Aerodynamics  1959-present 

Chairman,  Dept,  of  Mechanical  Engineering  1969-1972 

Visiting  Research  Officer  1982 


R.A.A.F. 

N.R.C.,  Canada 
Cambridge  University 
Miss.  State  University 
Univcrsite  Laval 
McGill  University 
McGill  University 
N.R.C. 


ACADEMIC  ACTIVITIES 


COURSES  TAUGHT 


305-33 1A-B 
305-533A 
305-610A 
305-612A 


Fluid  Mechanics 
Subsonic  Aerodynamics 
Fundamentals  of  Fluid  Mechanics 
Viscous  Flow  and  Boundary-Layer  Theory 


CONTRIBUTIONS  TO  UNIVERSITY,  FACULTY,  AND  DEPARTMENTAL 
GOVERNANCE. 


1961-1964 

1968-1972 

1973- 1975 

1974- present 

1975- 1979 

1977- 1979 

1978- 1979 

1979- 1982 

1980- present 


Secretary  of  Graduate  Faculty. 

Member  of  Committee  on  Research. 

Member,  Planning  Committee  of  Graduate  Faculty. 

Pro  Dean,  Graduate  Faculty. 

Member,  Tenure  Committee,  Faculty  of  Engineering. 
Member,  Nominating  Committee  for  Graduate  Faculty. 
Member,  McGill  Athletics  Review  Board. 

Member,  Senate. 

Member,  Nominating  Committee  of  Senate. 
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1989-present 


Member,  Council,  Faculty  Club:  Vice  President  1991-93; 

President  1993-95. 

Chairman,  Liaison  Committee  for  the  M.Eng.  in  Aerospace  Engineering. 


1989-present 


CONSULTING  AND  OTHER  PROFESSIONAL  ACTIVITIES 

Consultant,  Aircraft  and  aerodynamic  design  1959-1970 

Canadair  Ltd. 


Consultant, 

Expo  Corporation 

Wind  effects  on  buildings 

Consultant, 

Boeing  Scientific 
Research  Laboratories 

Jets  and  wall  jets 

Consultant, 

Dept,  of  Transport 

Licensing  of  light  aircraft 

Consultant, 

Asbestos  Corp. 

Wind  effects  on  the  roof  of  a 

German  building 

Consultant, 

Prau  &  Whitney 

Canada  Inc. 

Flow  and  heat  transfer  behind  impeller 
disks,  and  internal  aerodynamics,  jet 
flap  and  Coanda  jet  flap  cascades. 

Consultant, 

S.O.D.E.M.  Quebec 

Wind  turbines 

Consultant, 

Pulp  &  Paper 

Research  Institute 

Disk  refining  of  wood  pulp 

Consultant, 

Midland  Ross 

Analysis  of  web  stability  in  a  wood-pulp 
dryer 

1966 

1967 

1971 

1973 

1973-1991 

1975 

1975-1977 

1980-present 

1977-1980 


AWARDS  AND  HONOURS 

1961  Edward  Busk  Memorial  Prize  -  Royal  Aeronautical  Society. 

1964  Tour  of  British  Universities  and  Research  Establishments  under  the  auspices  of 

the  British  Council. 

1966  Invited  paper  General  Motors  Conference  on  ’Fluid  Mechanics  of  Internal  Flow’ . 

1967  N.R.C.  Senior  Travelling  Fellowship  -  Sabbatical  leave  in  Cambridge  University, 
Dept,  of  Theoretical  Physics  and  Applied  Mathematics. 

1969  Turnbull  Lecturer,  Canadian  Aeronautics  and  Space  Institute. 
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1975-1976  N.R.C.  Senior  Travelling  Fellowship  -  Sabbatical  leave  in  Cambridge  University 
Engineering  Dept. 


1975-1976  Visiting  Scholar,  Corpus  Christi  College,  Cambridge. 

1981  Invited  lecturer  CANCAM  ’81,  Moncton. 

1982  Visiting  Research  Officer,  National  Aeronautical  Establishment,  N.R.C. 

1987  Honorary  Chairman  CANCAM,  1987,  Edmonton. 


MEMBERSHIP  ON  NATIONAL  AND  INTERNATIONAL  BODIES 
1964  Fellow,  Royal  Aeronautical  Society. 


1964 

1965-present 

1961-64,  1967-69 
1978 

1967-1975 

1968 

1970- 1975 

1971- 1973 

1974-1976 

1984-1987 

1978- 1979 

1979- 1984 

1981 

1982 

1983 -present 

1983- 1985 

1984- 1987 

1985- 1987 
1987 

1993-1995 


Fellow,  Canadian  Aeronautics  and  Space  Institute. 

Member,  Order  of  Engineers,  Quebec. 

Member  and  Chairman,  N.R.C.  Associate  Committee  on  Aerodynamics. 

Member  and  latterly  Chairman,  D.R.B.  Advisory  Committee  for  Plasma 
and  Fluid  Dynamics. 

Chairman,  Technical  Programme  Committee,  C.A.S.I.  International 
Aerospace  Exposition. 

McGill  Member,  Canadian  Research  Management  Association. 

Member,  National  Committee  for  International  Union  of  Theoretical  and 
Applied  Mechanics. 

Member,  Council  Canadian  Aeronautics  and  Space  Institute. 

Coordinator  and  chief  reviewer  in  Fluid  Mechanics,  CANCAM  ’79. 
Member,  Aeronautics  Advisory  Board,  Dept,  of  Transport. 

Member,  A.S.M.E. 

Fellow,  Royal  Society  of  Canada. 

Member,  Academy  of  Science  of  the  Royal  Society  of  Canada. 

Member,  Advisory  Committee  on  Aerodynamics  for  the  N.A.E. 

Member,  Council  of  the  Royal  Society  of  Canada. 

Vice-President,  Academy  of  Science  of  the  Royal  Society  of  Canada. 

Founding  Fellow,  Canadian  Academy  of  Engineering. 

Director,  Division  of  Applied  Science  and  Engineering,  Academy  of 
Science  of  the  Royal  Society  of  Canada. 
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Curriculum  Vitae 
Theodore  Yao-tsu  Wu 

California  Institute  of  Technology,  Pasadena,  CA  91125  U.S.A. 
Telephone:  (818)  395-4230  Fax:  (818)795-9839 

e-mail:  tywu@romeo.caltech.edu 


Personal 

[Pll  Redacted] 

Education 

B. S.  (Aeronautics)  Chaio-Tung  University,  Shanghai,  China  1946 
M.S.  (Aeronautics)  Iowa  State  College,  Ames,  Iowa  1948 
Ph.D.  (Aeronautics)  California  Institute  of  Technology  1952 

Professional  Experience 

Faculty  member,  Division  of  Engineering  and  Applied  Science,  California  Institute 
of  Technology,  Pasadena,  CA,  U.S.A.,  1952- 
Research  Fellow,  1952-55; 

Assistant  Professor  of  Applied  Mechanics,  1955-57; 

Associate  Professor  of  Engineering  Science,  1957-61; 

Professor  of  Engineering  Science,  1961- 

Awards  and  Honors 

Guggenheim  Fellow,  1964 

Elected  member,  National  Academy  of  Engineering,  1982- 

Elected  member,  Academia  Sinica,  1984- 

Visiting  Professor,  Hamburg  Univ.,  West  Germany,  1964-65 

Achievement  Award,  Chinese  Engineers  and  Scientists  Assoc.  So.  CA,  1971 

C. S.I.R.O.  Fellow,  Australian  Universities  and  C.S.I.R.O.  Fellowship,  1976 
Adjunct  Professor,  Shanghai  Jiao  Tong  Univ.,  Shanghai,  China,  1979- 
Honorary  Professor,  Northwestern  Polytechnical  Univ.,  Xian,  China,  1979- 
Russell  Severance  Springer  Visiting  Professor,  UC,  Berkeley,  1980 

JSPS  Fellow,  Japan  Society  for  the  Promotion  of  Science  Fellowship,  1982 
Honorary  Professor,  Harbin  Shipbuilding  Eng.  Institute,  China,  1987- 
Honorary  Fellow,  Institute  of  Mechanics,  Academia  Sinica,  1988- 
Distinguished  Lecturer,  Hong  Kong  Univ.  of  Science  and  Technology,  1992 
Invited  Visiting  Professor,  Hong  Kong  Univ.,  Hong  Kong,  April  1992 
1993  CAFA  Achievement  Award  -  Chinese-American  Faculty  Assoc.,  So.  CA. 

1993  Fluid  Dynamics  Prize  Award  -  American  Physical  Society 

Reviewer  and  Referee  Service 

Journal  of  Fluid  Mechanics;  Physics  of  Fluids;  Mathematical  Reviews;  Publications 
of  ASME,  ASCE,  AIAA,  SNAME,  and  other  scientific  journals  in  addition  to  those 
journals  charged  with  editorship.  Reviewer  of  proposals  for  the  National  Science 
Foundation,  USARO(D)  and  other  agencies. 
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Publications 

Author  or  coauthor  of  over  100  professional  papers,  books  and  reports. 
Affiliations 

Sigma  Xi;  Phi  Tau  Phi;  Pi  Mu  Epsilon 

Member,  U.S.  National  Academy  of  Engineering,  Special  Fields  and  Interdisci¬ 
plinary  Engineering  Peer  Committee,  1992- 
Member,  Ges.  Angew.  Math.  u.  Mech.  (W.  Germany  &  Austria),  1965- 
Fellow,  American  Physical  Society,  1967- 

Vice-Chairman,  Division  of  Fluid  Dynamics,  American  Physical  Society,  1990-91; 

Member,  Executive  Committee,  Division  of  Fluid  Dynamics,  1991- 
Member,  Society  of  Naval  Architects  and  Marine  Engineers,  1972- 
Member,  Society  for  Industrial  and  Applied  Mathematics,  1987- 
Member,  Intern.  Towing  Tank  Conference,  Wave  Resistance  Committee,  1967-78; 
Chairman,  1977-78 

American  Towing  Tank  Conference,  1956-;  Chairman,  1972-74 
Hydrodynamics  H-5  and  H-8  Panels  of  Soc.  Nav.  Arch.  Mar.  Eng.,  1963 
Member,  AIAA,  Fluid  Mechanics  Committee,  1949-58 

National  Research  Council,  Committee  on  Recommendation  for  U.S.  Army  Basic 
Scientific  Research,  1975-78 

Program  Review  Counselor  for  the  College  of  Engineering,  University  of  Iowa, 
1976 

Promotion  Review  Committee,  Institute  of  Applied  Mechanics,  National  Taiwan 
University,  1990 

Consulting  Experience 

Aerojet  General,  Bendix  Corporation,  Chevron  Corporation, 

CONVAIR  General  Dynamics,  Douglas  Aircraft  Company,  General 
Electric,  Lockheed  Aircraft  Company,  North  American  Aviation, 

TRW  Systems  and  other  industrial  companies. 

Consultant  for  NRAC,  Dept,  of  the  Navy,  1982-83 

Consultant  for  David  Taylor  Research  Center,  Dept  of  the  Navy,  1985 

Editorship 

Co-editor,  Advances  in  Applied  Mechanics,  1982- 
Editorial  Board,  Advances  in  Applied  Mechanics,  1971- 
Editorial  Committee,  Journal  of  Ship  Research,  1972- 
Editorial  Board,  The  Physics  of  Fluids,  1979-82 
Editorial  Board,  Wave  Motion,  1978- 

Editorial  Committee,  Annual  Review  of  Fluid  Mechanics,  1980-84;  1691- 
Editorial  Committee,  Applied  Mathematics  and  Mechanics,  1979- 
Editorial  Committee,  Newton  -  Graphic  Science  Magazine,  1986- 
Editorial  Committee,  Advances  in  Hydrodynamics,  1988- 
Advisory  Board,  Journal  of  Acta  Mechanics  Sinica,  1990- 
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John  L.  Lumley 

B.A.  -  Harvard  College,  1952,  Engineering  Sciences  and  Applied  Physics 

M.S.E.  -  The  Johns  Hopkins  University,  1954,  Mechanical  Engineering 

Ph.D.  -  The  Johns  Hopkins  University,  1957,  Aeronautics.  Particular  emphasis  on 

turbulence  and  stochastic  processes. 

Haute  Distinction  Honoris  Causa  -  Ecole  Central  de  Lyon,  1987. 

Fellow,  American  Academy  of  Arts  and  Sciences. 

Member,  National  Academy  of  Engineering. 

Timoshenko  Medal,  American  Society  of  Mechanical  Engineers,  1993 
ALAA  Fluid  and  Plasmadynamics  Award,  1982. 

APS  Fluid  Dynamics  Prize,  1990. 


Post-doctoral  fellow,  J.H.U.,  1957-1959.  The  Pennsylvania  State  University  at 
University  Park,  1957-1977;  F*rofessor  of  Aerospace  Engineering,  1963-74;  Evan  Pugh 
FWessor  of  Aerospace  Engineering,  1974-1977.  Cornell  University  (Sibley  School  of 
Mechanical  and  Aerospace  Engineering),  Willis  H.  Carrier  F*rofessor  of  Engineering, 
1977-present.  Fulbright  Senior  Lecturer,  Belgium;  Guggenheim  Fellow,  France,  1973-74; 
Author,  Stochastic  Tools  in  Turbulence  (1970);  co-author,  Structure  of  Atmospheric 
Turbulence  (with  H.A.  Panofsky;  1964);  A  First  Course  in  Turbulence  (with  H. 

Tennekes;  1971).  Principal,  Deformation  of  Continuous  Media  (1963),  Eulerian  vs. 
Lagrangian  Frames  in  Fluid  Mechanics  (1968)  (Motion  Picture  Films).  Technical  Editor, 
Statistical  Fluid  Mechanics  (by  A.S.  Monin  and  A.M.  Yaglom)  (1971,  1975);  Variability  of 
the  Oceans  (by  A.S.  Monin,  V.M.  Kamenkovich  and  V.G.  Kort)  (1977).  Member,  Johns 
Hopkins  Society  of  Scholars,  Sigma  Xi,  New  York  Academy  of  Sciences,  Society  for 
Natural  Philosophy,  American  Association  for  the  Advancement  of  Science;  Associate 
Fellow,  American  Institute  of  Aeronautics  and  Astronautics;  Fellow,  American  Physical 
Society  (Division  of  Fluid  Dynamics),  American  Academy  of  Mechanics. 

Research  papers  (Last  Five  Years) 


124.  Shih,  T.-H.  &  Lumley,  J.  L.  1986.  Second  order  modeling  of  particle  dispersion 
in  a  turbulent  flow.  J.  Fluid  Mech.  163:  349-363 

126.  Shih,  T.-H.  &  Lumley,  J.  L.  1986.  A  general  form  of  the  equations  for  turbulent 
flows  with  a  passive  scalar.  Mathematical  Modeling.  7:  353-370. 

128.  Shih,  T.-H.,  Lumley,  J.  L.  &  Janicka,  J.  1987.  Second  order  modeling  of  a 
variable  density  mixing  layer.  J.  Fluid  Mech.  1 80:  93- 1 1 6. 

133.  Lumley,  J.L.  1986.  Evolution  of  a  non-self-preserving  thermal  mixing  layer. 
Phys.  Fluids  29  (12):  3976-3981. 

134.  Leibovich,  S.  &  Lumley,  J.L.  1986.  Complex  fluid  motions:  Models  and 
metaphors.  Cornell  Engineering  Quarterly  20:  27-35. 

135.  Lumley,  J.L.  1987.  Turbulence  modeling.  In  Proceedings  of  the  Tenth  U.S. 
National  Congress  of  Applied  Mechanics ,  ed.  J.P.  Laub.  pp.  33-39.  New  York: 
ASME. 
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136.  Aubry,  N.,  Holmes,  P.,  Lumley,  J.L.  and  Stone,  E.  1987.  The  dynamics  of 
coherent  structures  in  the  wall  region  of  a  turbulent  boundary  layer.  J.  Fluid  Mech. 
In  press. 

137.  Aubry,  N.,  Holmes,  P.,  Lumley,  J.L.  and  Stone,  E.  1987.  Models  for  coherent 
structures  in  the  wall  layer.  In:  Advances  in  Turbulence,  eds.  G.  Comte-Bellot,  J. 
Mathieu,  pp.  346-356.  Berlin/Heidelberg:  Springer. 

138.  Turbulent  Shear  Flows  5,  eds.  F.  Durst,  B.E.  Launder,  J.L.  Lumley.  F.W. 
Schmidt,  J.H.  Whitelaw.  Heidelberg:  Springer. 

139.  Sarkar,  S.  &  Lumley,  J.L.  1987.  Angular  dispersion  of  material  lines  in  isotropic 
turbulence.  Phys.  Fluids.  30  (5):  1269-1271. 

140.  Lumley,  J.  L.  1987  The  state  of  turbulence  research.  In  Advances  in  Turbulence, 
eds.  W.  K.  George  &  R.  E.  A.  Arndt,  pp.  1-10.  New  York:  Hemisphere. 

141.  Aubrey,  N.  Holmes,  P.H.,  Lumley,  J.L.  and  Stone,  E.  1990.  The  behavior  of 
coherent  structures  in  the  wall  region  by  dynamical  systems  theory.  In  Near-Wall 
Turbulence  eds.  S.J.  Kline  &  N.  H.  Afgan,  pp.  672-691.  Washington,  DC: 
Hemisphere. 

142.  Aubry,  N.,  Holmes,  P.,  Lumley,  J.L.  and  Stone,  E.  1989.  A  simple  model  for 
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Professor  S.  C.  R.  Dennis  and  concerned  the  development  of  calculation  methods  for  the 
viscous  steady  and  unsteady  flow  past  spheres. 

Upon  completion  of  his  Ph.D.  studies  in  1971,  Dr.  Walker  won  a  National  Research  > 

Council  of  Canada  Postdoctoral  Fellowship  (the  equivalent  of  a  NATO  Postdoctoral 
Fellowship)  to  work  with  the  late  Professor  K.  Stewartson,  F.R.S.  at  University  College 
London  in  the  University  of  London,  England.  During  the  period  from  1971  -  1973,  he 
carried  out  research  related  to  compressible  unsteady  boundary-layer  flows  in  shock 
tubes,  flow  separation  in  rotating  and  magnetohydrodynamic  flows,  and  heat  transfer 
processes  for  low  Reynolds  number  flows  past  spheres. 

In  1973,  he  joined  the  School  of  Mechanical  Engineering  at  Purdue  University  as  an 
Assistant  Professor.  During  this  period,  he  carried  out  research  in  the  modeling  of 
turbulence  and  heat  transfer  processes  in  boundary  layers,  in  unsteady  laminar  flows 
and  in  the  evaluation  of  the  heat  and  mass  transfer  processes  that  take  place  in  rapidly- 
rotating  gas  centrifuges. 

Dr.  Walker  joined  Lehigh  University  in  1978  as  an  Associate  Professor  and  was 
subsequently  promoted  to  Professor  of  Mechanical  Engineering  and  Mechanics  in  1983. 

His  research  at  Lehigh  has  been  associated  with  turbulence  modeling,  heat  and  mass 
transfer  processes  in  turbulent  boundary  layers  and  in  rapidly  rotating  gas  centrifuges, 
unsteady  viscous  boundary-layer  flows,  numerical  methods,  the  prediction  of 
electroplating  rates  on  metal  surfaces,  and  hypersonic  boundary-layer  separation.  He  is 
nationally  recognized  for  his  research  in  unsteady  flows  and  turbulence  and  over  the 
past  several  years  has  given  several  invited  presentations  at  various  international 
meetings,  workshops,  and  national  scientific  meetings. 

Over  the  years,  Dr.  Walker  has  been  the  principal  or  co-principal  investigator  on  thirty- 
one  research  grants  and  contracts  totaling  approximately  $6,500,000.  These  include 
several  grants  and  contracts  from  NASA,  AFOSR,  United  Technologies  and  AMP 
Incorporated.  Dr.  Walker  has  served  as  the  Project  Director  of  a  combined 
Industry/University  NSF  grant  with  United  Technologies  Research  Center  to  study 
heat  transfer  effects  on  gas  turbine  blades  and  as  a  Principal  Investigator  on  a  recent 
AFOSR  URI  entitled  “Computational  and  Analytical  Methods  in  Nonlinear  Fluid 
Mechanics”.  He  is  the  author  of  over  eighty  refereed  technical  journal  publications, 
conference  proceedings  and  other  publications  in  the  areas  of  fluid  mechanics,  heat  and 
mass  transfer,  and  numerical  methods.  Twenty-nine  graduate  students  (fourteen 
Ph.D.’s)  have  graduated  with  Dr.  Walker  as  their  major  professor.  He  currently  advises 
six  Ph.D.  candidates  and  four  Master’s  students. 

He  is  a  member  of  AIAA,  the  American  Physical  Society,  Sigma  Xi  and  ASME  and  was 
selected  to  the  AIAA  Fluid  Dynamics  Technical  Committee  in  1989.  He  is  a  Fellow  of 
the  American  Physical  Society  and  an  Associate  Fellow  of  AIAA.  He  was  elected 
Chairman  of  the  Fluid  Dynamics  Technical  Committee  in  1992  for  the  term  1993-1996. 

Dr.  Walker  is  currently  Chairman  of  the  Thermo-Fluids  Division  in  the  Department  of 
Mechanical  Engineering  and  Mechanics  at  Lehigh  University. 


A-16 


Recent  Representative  Publications  (Total  Publications  -  85) 

1.  Doligalski,  T.  L.,  Smith,  C.  R.  and  Walker,  J.  D.  A.  1994  “Vortex  interactions 
with  walls”,  Ann.  Rev.  Fluid  Mechs.  26,  pp.  578-616. 

2.  Kerimbekov,  R.  M.,  Ruban,  A.  I.  and  Walker,  J.  D.  A.  1994  “Hypersonic 
boundary-layer  separation  on  a  cold  wall”,  J.  Fluid  Mech.  (in  press). 

3.  Degani,  A.  T.,  Smith,  F.  T.  and  Walker,  J.  D.  A.,  “Asymptotic  structure  of  three- 
dimensional  turbulent  boundary  layers”,  J.  Fluid  Mech.  250,  pp.  43-68,  1993. 

4.  Smith,  C.  R.,  Walker,  J.  D.  A.,  Haidari,  A.  H.  &  Sobrun,  U.,  “On  the  dynamics  of 
near-wall  turbulence”,  Phil.  Trans.  Roy.  Soc.  Lond.  A  336,  pp.  131-175,  1991. 

5.  He,  J.,  Kazakia,  J.  Y.,  Ruban,  A.  I.  &  Walker,  J.  D.  A.,  “An  algebraic  model  for 
dissipation  in  supersonic  turbulent  boundary  layers”,  AIAA  Paper  92-0311. 

6.  Degani,  A.  T.  &;  Walker,  J.  D.  A.,  “Computation  of  three-dimensional  turbulent 
boundary  layers  using  embedded  functions”,  AIAA  Paper  91-0440,  submitted  to 
AIAA  J. 

7.  Degani,  A.  T.,  Walker,  J.  D.  A.,  Power,  G.  &  Ersoy,  S.  “On  the  application  of 
algebraic  turbulence  models  to  high  Mach  number  flows”,  AIAA-91-0616. 

8.  Degani,  A.  T.,  Smith,  F.  T.  fc  Walker,  J.  D.  A.,  “The  three-dimensional  turbulent 
boundary  layer  near  a  plane  of  symmetry” ,  J.  Fluid  Mech.  234,  pp.  329-360,  1991. 

9.  Peridier,  V.  J.,  Smith,  F.  T.  &  Walker,  J.  D.  A.,  “Vortex-induced  boundary-layer 
separation,  Part  I:  The  limit  problem  Re  -»  oo”,  J.  Fluid  Mech .”  232.  pp.  99-131, 
1991. 

10.  Peridier,  V.  J.,  Smith,  F.  T.  &  Walker,  J.  D.  A.,  “Vortex  induced  boundary-layer 
separation.  Part  II  -  Unsteady  interaction  boundary-layer  theory”,  J.  Fluid  Mech. 
232.  pp.  133-165,  1991. 

11.  Hoyle,  J.  M.,  Smith,  F.  T.  &  Walker,  J.  D.  A.,  “On  sublayer  eruption  and  vortex 
formation”,  Comp.  Phys.  Commun.  6&>  pp.  151-157,  1991. 

12.  Walker,  J.  D.  A.,  Ece,  M.  C.  &  Werle,  M.  J.,  “An  embedded  function  method  for 
turbulent  flow  prediction”,  AIAA  Paper  87-1464,  19th  Fluid  Dynamics,  Plasma 
Dynamics  and  Lasers  Conference,  Honolulu,  Hawaii,  June  1987;  AIAA  J.  22,  pp. 
1810-1818,  November,  1991. 

13.  Hon,  T.-L.  &  Walker,  J.  D.A.,  “Evolution  of  hairpin  vortices  in  a  shear  flow”, 
NASA  Technical  Memorandum  100858,  ICOMP-88-9,  NASA  Lewis  Research 
Center,  Cleveland,  Ohio,  1988;  Comp.  &  Fluids  20,  pp.  343-358,  1991. 

14.  Walker,  J.  D.  A.,  “Models  based  on  dynamical  features  of  the  wall  layer”,  Appl. 
Mech.  Rev.  43,  No.  5,  Part  2,  pp.  5232-5239,  1991. 

15.  He,  J.,  Kazakia,  J.  Y.  &  Walker,  J.  D.  A.,  “Embedded  function  methods  for 
supersonic  turbulent  boundary  layers”,  AIAA  Paper  90-0306,  January  1990. 


A-17/A-18 
Reverse  Blank 


Curriculum  Vitae:  KATEPALLI  R.  SREENIVASAN 
Mason  Laboratory,  Yale  University,  New  Haven,  CT  06520-2150 

Tel:  (203)  432-4345,  Pax:  (203)  432-1054 
c-mail :  krs@kolmogorov.eng.yale.edu 


[Pll  Redacted] 


Academic  Degress:  /  A  „  VT 

B.E.  (Mech.  Engg.)  Bangalore  University  (I  Rank),  1908;  M.E.  (Aero.  Engg.)  Indian 
Institute  of  Science  (1  Rank),  1970;  Ph.D.  (Aero.  Engg.)  Jndiun  Institute  of  Science 
(gold  medal),  1975;  M.A.  (Privatim)  Yale  University,  1986 

CuirenlMiioa  ...  .  , 

Harold  W.  Cheel  Professor  ol  Mechanical  Engineering  (concuiTcnt  uppomunenls  ui  the 
Departments  of  Physics  and  Applied  Physics) 

Professional  Society  Membership: 

American  Academy  of  Mechanics;  American  Institute  of  Astronautics  and  Aeronautics; 
American  Mathematical  Society;  American  Physical  Society;  American  Society  of 
Mechanical  Engineers;  Society  for  Natural  Philosophy;  Sigma  Xi 

Awards  and  Honors: 

Sir  M.  Visvesvaraya  Memorial  Prize,  1968;  Indumati  Lalbai  Memorial  Gold  Medal, 
1968;  Institution  of  Engineers  Prize  and  Honorary  Graduate  Membership,  1968;  J.  R. 
D.  Tata  Fellow,  1973;  P.  S.  Narayana  Gold  Medal  for  (lie  best  Ph.D.  thesis  in 
Mechanical  Sciences,  1975,  Humboldt  Fellow,  1983;  Fellow,  American  Physical 
Society,  1985;  Elected  to  the  Connecticut  Academy  of  Science  and  Engineering,  1988; 
Guggenheim  Fellow,  1989;  Induction  to  Society  of  Scholars,  Johns  Hopkins 
University,  1991 ;  Distinguished  Alumni  Award  of  the  Aerospace  Department,  Indian 
Institute  of  Science,  1992;  Fellow,  American  Society  of  Mechanical  Engineers,  1993. 

A  summary  of  professional  and  administrative  activities: 

Publications  relating  to  turbulent  flows,  sonic  booms,  nuclealion  of  droplets  in 
condensation,  chaos,  fractals,  and  so  forth,  about  100. 

At  various  times,  associate  editor  or  a  member  of  the  editorial  board  of;  Journal  of 
Applied  Mechanics  (ASME);  Proceedings  of  the  Indian  Academy  of  Science,  scries  C 
(Sadhana);  American  Scientist;  Journal  of  Nonlinear  Science;  Springer  Book  Series  on 
Interdisciplinary  Applied  Mathematics;  Journal  of  Theoretical  and  Computational  Fluid 
Dynamics;  Physical  Review  Letters 

Delivered  numerous  invited  lectures  (including  mini-courses)  in  Universities,  Institutes, 
and  Government  Laboratories  in  the  U.S.,  Germany,  Great  Britain,  France,  Australia, 
India,  Poland,  Norway,  the  Netherlands,  Italy,  etc.  Among  them:  Emerging  Scholar 
lecture,  University  of  Notre  Dame;  the  First  Stanley  Corrsin  Memorial  Lecture,  Johns 
Hopkins  University;  Phillips  Lecture,  Havcrford  College;  Sabita  Choudhury  Memorial 
Lecture,  Indian  Institute  of  Science 

Invited  speaker  in  many  international  meetings  and  workshops  on  fluid  dynamics  and 
physics,  turbulence,  nonlinear  dynamical  and  fractals 
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Organized  several  scientific  workshops,  meetings  and  discussion  sessions  in  meetings 
sucji  as  IUTAM,  SIAM,  APS,  3rd  World  Conference  on  Experimental  Heat  Transfer, 
Fluid  Mechanics  and  Thermodynamics,  etc. 

At  various  limes:  Representative  of  die  Division  of  Fluid  Dynamics  of  the  American 
Physical  Society  for  the  Committee  on  the  International  Freedom  of  Scientists;  member 
of  the  Executive  Committee  of  the  Division;  member  of  the  Fluid  Dynamics  Prize 
Committee;  member  of  the  Frenkiel  Award  Committee;  Chairman  (und  the  Past 
Chairman)  of  the  Fluid  Dynamics  Division  of  the  American  Physical  Society;  member 
of  the  publications  committee  of  the  APS  (DFD);  member  of  the  Fluids  Engineering 
Division  of  the  American  Society  of  Mechanical  Engineers,  etc.. 

Participated  in  collaborative  atmospheric  and  oceanic  measurement  expeditions  in 
Australia,  as  well  as  in  field  experiments  on  sonic  booms  in  India 

At  Yale:  as  Chairman  of  die  Mechanical  Engineering  Department,  for  five  and  one-half 
years  generally  responsible  for  its  administration;  specifically,  for  appointments  and 
promotions  and  for  renewing  the  ABET-acereditation  of  the  Mechanical  Engineering 
Curriculum;  as  acting  Chairman  of  the  Council  of  Engineering  (Applied  Physics, 
Chemical  Engineering,  Electrical  Engineering  and  Mechanical  Engineering)  for  half-a* 
year,  responsible  for  the  management  of  the  Council  Administration;  member  of  the 
Advisory  Committee  for  Physical  Sciences  and  Engineering 

Miscellaneous: 

Mentioned  in  "American  Men  and  Women  of  Science",  "Who  is  Who  in  the  East," 
"Who  is  Who  in  the  World",  etc. 

Visiting  positions  held: 

Visiting  Professor,  Aerospace  Department,  University  of  Tucson,  Arizona;  Visiting 
Scientist,  Center  for  Atmospheric  Sciences,  Indian  Institute  of  Science;  Visiting 
Scientist,  DFVLR,  Goettingen,  West  Germany;  Visiting  Professor  of  Aeronautics, 
California  Institute  of  Technology;  Visiting  Professor  of  Physics,  Rockefeller 
University;  Visiting  Professor,  Jawaharlal  Nehru  Center  for  Advanced  Scientific 
research.  Bangalore,  India. 

Graduate  students  .supervised.,  and  post-doctoral  fellows  and  visitors: 

Fourteen  Ph.D.'s  completed;  about  a  dozen  post-doctoral  fellows  and  visitors  to-date 
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02244  (routing  copy) 

0251  (Conforti) 

0261 
0262  (2) 

10  (Nadolink) 

10B  (Griffin) 

102  (Lima) 

103  (Dickinson) 

104 
20 

213  (Schloemer,  Snarski,  Payne) 

214  (Bakewell,  Keith,  Hurdis,  Chen,  Abraham) 

2211  (Graham,  Nguyen) 

23132  (Keefer) 

38 

382 

383 

411  (Kalinowski,  Wagner) 

60 

72 

80  (Dence,  Myers) 

801  (Meng) 

814 

82  (White) 

8211  (Aaron) 

8212  (Shonting) 

8222  (Nugent) 

8213  (Impagliazzo) 

8214  (Kemp) 

8231  (Khoury) 

8232  (Hendricks,  Hrubes) 

8233  (Bandyopadhyay  (30),  Beauchamp,  Dick,  Grant,  Henoch,  Huyer, 

Philips,  Rice,  Uhlman,  Varghese) 

8234  (Corriveau,  Harari) 

83 

832 

8322  (Jordan) 

833 


Total: 
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